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A QUANTUM MECHANICAL DESCRIPTION OF 
COLLECTIVE MOTIONS! 


By R. SKINNER? 


ABSTRACT 


A rigorous and general formalism is presented for the quantum mechanical 
description of the collective motions of a system containing a finite number of 
articles. The dynamical equations are expressed in terms of collective and 
‘internal” coordinates. The internal coordinates are symmetrical in the particles, 
and correspond to the laboratory coordinates with the collective motions 
subtracted off. The formalism is applied to center of mass motion. 


INTRODUCTION 


The recent investigations by Bohr and Mottelson (1953), Hill and Wheeler 
(1953), and others, of the combined liquid drop — individual particle model 
of the nucleus have led to an increased interest in exact methods for obtaining 
the equations for collective motions of systems having a finite number of 
particles. The collective equations of motion may be derived by a point 


transformation from the 3A particle coordinates (x;, y;, 2;) to the N collective 
coordinates 6, and (3A —JN) ‘‘internal’’ coordinates. These internal coordinates 
are equivalent to (3A —JN) of the 3A particle coordinates. This straightforward 
procedure gives the dependence of physical variables, such as the energy, on 
the collective coordinates and the derivatives with respect to these coordinates. 
However, it destroys the symmetry which exists when these quantities are 
expressed in terms of all the individual particle coordinates. This asymmetry 
increases the complexity of the calculations. Moreover, physical variables 
which depend only on the (3A — J) particle coordinates occur as the coefficients 
of the collective coordinates and their derivatives in the dynamical equations 
for the collective motions. It is difficult to see the dependence of these physical 
variables on the N missing particle coordinates. 

Lipkin, de Shalit, and Talmi (1955) have introduced a formalism in which 
symmetry is preserved by introducing N superfluous coordinates, and also N 
superfluous degrees of freedom. The present formalism preserves the symmetry 
by introducing N superfluous coordinates, and yet corresponds exactly to the 
original system by not changing the number of degrees of freedom. Others 
(e.g., Bohr 1954; Tolhoek 1955; and Tomonaga 1955) have considered the 


i1Manuscript received in original form December 21, 1955, and, as revised, May 20, 1956. 
Contribution from the Physics Department, University of Saskatchewan, Saskatoon, 


Saskatchewan. 
*Holder of a National Research Council of Canada University Postdoctorate Fellowship. 
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problem, but have not presented an exact formalism using internal coordinates 
symmetrical in the particles. The formal treatment presented here is rigorous 
and provides a general method for describing the collective and internal 
motions of many-particle systems. The treatment may be used to show the 
extent of the effect of approximations that are made in less exact treatments. 

Let the collective motions under consideration be described by the coordi- 
nates 6;,..., Oy, where N is the number of degrees of freedom of these motions. 
That part of the motion which is not included in the collective motions will be 
described by the internal coordinates. These are position vectors for each 
particle, 01, ..., @4, in a coordinate system equivalent to the laboratory system 
with the collective motions subtracted off. These 3A internal coordinates 
correspond to (3A —J) of the degrees of freedom and so must be subject to V 


constraints 
¢n(o:) — 0, a= iL. cee, N. 


A system of particles with laboratory coordinates 9; satisfying these constraints 
cannot undergo the collective motions. For example, the relation 


dL mio: = 0, 


satisfied by the internal coordinates considered in center of mass motions, is 
equivalent to the restriction on a system of particles with laboratory co- 
ordinates 9; that the center of mass remain fixed at the origin. 

A symmetric (1, 1) transformation connecting the 3A-dimensional subspace 
of possible motions in the (34+ )-dimensional space of the collective and 
internal coordinates with the 3A-dimensional space of the laboratory coordi- 
nates is obtained as follows. We assume that the transformation 


x; = X:(0;, On); i= A. sessile 


is known. Then the laboratory configuration space is extended to (83A+JN) 
dimensions through the introduction of NV new coordinates 


Xn = n(0;), n= Aes 


We restrict ourselves to the original 3A-dimensional space by defining each x, 

to be zero. The 3A-dimensional momentum space is similarly extended through 

the definition of NV new ‘‘momenta’”’ p, which are also defined to be zero. The 

Poisson brackets of the new quantities x, and p, with any other expression 

must be zero because, when we quantize the dynamical system, we must have 
[xn, f] = th(Xnf—fxn) = 0 = [pny f] 

for any f. 

The constraint functions ¢, are not arbitrary. Their functional form depends 
on the chosen collective coordinates 6, since it is necessary that the Jacobian 
of the transformation 

Xi, dn > Oi, On 


be non-zero but finite. Furthermore, the collective coordinates are useful only 
if matrix elements of operators like O;(9;)02(0,) between states with separable 
state functions x(0;)¢(6,) are of the form 


{fxa*Oixegi(o:)d%01..- }{fba*Ordoge(On)d01...}. 
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Therefore it is desirable that the Jacobian be separable into the product of a 
factor depending on the internal coordinates only and one depending on the 
collective coordinates only. 

The formalism is applied to the treatment of the simplest of collective 
motions—center of mass motion. 

The details of the calculations are simplified by assigning each coordinate 
in configuration space a single subscript 7 = 1, 2,..., 4f = 3A. The sub- 
scripts may be considered to refer to particles if the coordinates are under- 
stood to be vectors and the derivatives to be gradient operators. The elements 
of the matrix in Appendix A will be 3X3 matrices if the number of collective 
coordinates is divisible by three. 


CLASSICAL THEORY 


Let x; (¢ = 1, 2,..., M) be the laboratory coordinates of the particles in a 
dynamical system characterized by the Lagrangian 


(1) Ee wBlde ns ix Mens Be 


where © = d/dt. The collective motions are described by the collective coordi- 
nates p, (n = M+1,..., M+N).* The internal coordinates p; coincide with 
the laboratory coordinates when the collective coordinates are all zero. These 
internal coordinates are prevented from undergoing the collective motions by 
the NV constraints 


(2) on(p1, coos pm) = 0. 


The transformation equations connecting the laboratory coordinates with the 
collective and internal coordinates are 


(3) Xi = tulpn +++» Ox) 
The asymmetry in the number of coordinates in the two systems is removed by 
introducing N new “‘laboratory coordinates’”’ x, through the definitions 
(4) Xn = on(pr, ae ere yg Pa) _ 0. 

The dynamical system can be specified, in the laboratory coordinate system, 
by the Hamiltonian H, a function of the coordinates x; and their conjugate 


momenta p;, together with the Poisson bracket relations between these 
dynamical variables : 


(5) pi = OL /dx;, 
(6) H = px;-L, 
(7) [x:, x;] = 0, [x pi] = 653, [pi p3) = 0, 


where 6;; is the Kronecker delta. 
The transformed system is characterized by a Lagrangian &, 


(8) GS ae Bip, «+. 5 Rt Bae se Bes 
(8’) L(p1, sey px) = L[x1(1, eoeey px), a | Xu(p1, sees px)] ’ 

*The letters... , i, 7, k, / are restricted to the values 1,..., M throughout; m, n,... to the 
values M+1,..., M+N; and a, 8,... to the values 1,..., M+N = X. Summation is 


understood over repeated indices over the range of values specified by those indices. 
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but may also be specified by a Hamiltonian § and the Poisson bracket relations 
between the dynamical variables. Momenta canonically conjugate to the p; 
cannot be introduced because of the equations (2). We use as dynamical varia- 
bles the coordinates p, and variables that we shall call ‘‘momenta belonging to 
the coordinates’’, 7,. These momenta are defined through the relations 


(9) Ta = OL/Opa; 
and are related to the p; through the equations 
(10) Ta = O%;/Opa . Pi = %;/Apa » Pi. 
We remove the asymmetry in the number of momenta in the two systems 


by introducing NV new “laboratory momenta belonging to the x,’’, pn, which 
are defined to be zero: 


(11) O = Pa = xn(p1,.--» Px, T1y-++» Wx)» 
The functions x, are determined by rewriting (10) as 

(12) Ta = OXg/Opa - Pp 

This equation is equivalent to the relation 

{13) bs = 9pa/dxg . Ta, 


where || dp./dxg || is the matrix inverse to ||0x2/dp,|| . It may be calculated by 
solving equations (3) and (4) for the p, and then differentiating the solutions 
or, rather than solving these (possibly transcendental) equations, by inverting 
the X by X matrix ||0x./dp, ||. The latter procedure may be reduced to the 
inversion of an NXWN matrix if each transformation function x;(p1,..., px) 
depends only on the corresponding p; and the collective coordinates. This is 
shown in Appendix A. (If x; and p; are vectors, then it will be necessary that 
N be divisible by three and that the elements 0x./dpg be considered 3X3 
matrices for the inversion method of Appendix A.) Equations (11) and (13) 
show that 
(14) Va = Opa/OXn Ka = 0: 
This equation also results from the fact that the Lagrangian is independent 
of 2,. 

The new variables x, and p, have zero Poisson brackets with every dy- 
namical variable : 
(15) [xa, Xn] = [Xa Pn] = [xn, Pal = [Pn Pa] = O. 
The Poisson brackets of the collective and internal dynamical variables are 
calculated from the bracket expressions (7) and (15) in Appendix B: 

[Pas pa] = 0, [Pas 7] = (Opa/ 9x5) (dx ;/Op~), 
(16) [tay ™s] = (0x;/Opp)(O{ Ox ;/Opa} /Ax;)(Ppy/dx;) my 
— (0% 5/Opa)(O{ Ax ;/Apg} /Ax;:)(Ap_/dx;) my. 


The bracket expressions could have been deduced by using the brackets of 
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Dirac (1950). This procedure involves a knowledge of the functions x, and the 
inversion of a2NX2N matrix. 

In the quantized system, we shall want coordinate and momentum operators 
which are separable into two sets, one containing operators for the collective 
motions only and the other operators for the internal motions only. Further- 
more, it is necessary that (a) each operator of the first set have zero Poisson 
bracket with each operator of the second set and (b) every operator xj, p; 
having physical meaning in the laboratory system be a function only of the 
members of these two sets in the collective and internal coordinate systems. 
The first condition is necessary so that separable wave functions 


W = x(p1,.--, par) O(paret, - ++ > Pmtn) 


can certainly exist (Dirac 1947, §20). This condition is not satisfied by the sets 
{pi, w;} and {pm, mn} since, for example, the Poisson bracket 


len, 7j] = (Apn/Ax;)(Ox;/dp;) = —(Apn/IXm)(dXm/dp;) 
is not zero in general. However, the M linear combinations of 7, 
(17) Oo; = Ti —Ai;7;, 


have zero Poisson bracket with p, if the a;; satisfy the relations 


(18) OXn/Opi—i;(0Xn/9p;) = 0. 
‘Since we want the o; to depend on the internal coordinates and momenta only, 
we impose the condition that a;; be independent of py41,..., pain. Now, 
ai;' = (0X%m/0p;:)(Ap;/dxXm) is a solution of equation (18), for 
(19) (0p; /OXm) (dxn/9p;) = (Opa/AXm) (0Xn/Opa) 

a (O54: 
This relation follows since x, is independent of the variables pyy41,... , pwn 


Although the (dx,/dp;) are functions of the internal coordinates only, in 
general the (0p;/dxm), and therefore the a;;’, depend also on the collective 
coordinates. We note that equation (19) is true for all values of the collective 


coordinates, and in particular for py4i1 = ... = px = 0. Therefore, an accept- 
able solution for a;; is 
(20) iz = (a:;')o = (AXm/dp;)(0p;/AxXm)o, 
where 
( )o - ( ) a or 


The corresponding momentum operator is 
(21) Oi = 7i—(OXm/Op;) (Op; /OXm) 0}. 


The sets of variables {p;, o;} , {pn, mm} satisfy condition (a), for, as shown in 
Appendix B, the Poisson brackets of these variables are 


(22) [Ons Pm] = 0, [Ons Tm] = Sam) [n, Tm] = 0, 


[Pns pi] _ [ tn, pi] = [Pns a4] — [mn a1] = 0, 
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[oi es] = 9, — [pu o 5] = 51y— (OXm/Op;) (Op:/AXm)o, 
[o 1, 0)] = 0/Opi{ (OXn/Op;) (Opx/IXn)o} o% 
(23) ( “" —8/8p;{ (Ax%n/Ap:) (Ope /AXn)o} ox 
| — (%m/p1)(dpe/A%m)o 8/Ape{ (8%n/ ps) (Bpn/%n)o}on 
+ (OXm/Op;) (Ope /IXm)o O/Apx{ (OXn/ Ap 1) (Apn/AXn)o} on 


The sets {p:, o;} , {pn, tm} also satisfy condition (b) since the laboratory 
momenta, p;, are given by 


(24) pi = (0p;/0x;)0;+(Opn/dx:) tn. 
This equation is derived from equations (13) and (21) and the relation 
(Ap;/Ox;)(OXm/Op;) = (Opa/Ox;)(OXm/Opa) = O. 


The momenta a; are defined by equation (24) and a relation similar to (11). 
The functions x, of (11) cannot be expressed in terms of o; and 7,. The 
equations x, = 0 are replaced by the equations 
(25) (0p;/OXn)oo; = 0. 


The left-hand side of (25) could be used as the definition of p, rather than 
equation (11). Equations (24) and (25) define 7, and o;, so that (25) is satisfied 
identically by the internal variables o;, p;. Note that equations (25), like (2), 
do not contain the collective coordinates or momenta. 

The bracket expressions (23) may be simplified by using the fact that 
x; = p; when the collective coordinates are zero. The functions x;(p;, pn) may 
be expanded in a Taylor’s series with respect to pasy1,..-, Px! 


x= Pit (Oxx/Opn)o Pn + sae 


where the omitted terms contain a factor of p, at least squared. Since the pp 
are independent of p;, we have 


(Oxx/Op;)o = 4x;. 


Because 
(0X m/Op;) (0p;/OXm)o = bi;— (0xx/Op;)o (0p:/Axx)o, 


equations (23) may be written 
[ops] = 9, — [os 75] = (0p: /0X5)o, 
(23’) (71, 75] = (px /Ox;)0 9/Apx{ (Apn/Ox:)o}on 
— (Opx/ Ax :)0 O/Opx{ (Opn /Ox;)o} on. 


The Hamiltonian § in the new system can be derived from the definition 
(cf. Dirac 1950) 


D = patTa—% 


or, equivalently, 
(26) D = Pio + Pntn—%. 
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The velocities may be eliminated from (26) by solving the momenta-defining 
equations 

(0p; /0x;)0;+(Opn/dx;)an = AL/0%; 
for the combinations (0x;/0pa)pa. Alternatively, (26) may be rewritten as 

© = X,[(0p;/dx;)0;+(Apn/dx;)rn] —L 

= %;p;—L. 

The velocities can be eliminated from this simply by substituting for the x;, p; 
in the solutions for x; of equations (5). Therefore 
(27) D(p1,..., px, 01,..., 8x) = Alxi(pi,..., px),---> Pulor---, Fx). 


The equations of motion 
(28) a = [a, 9] 
follow from relation (27) and the definition of the Poisson brackets of the 
collective and internal variables. 
QUANTIZED THEORY 

The classical dynamical system is quantized by replacing the dynamical 
variables with Hermitian operators and the Poisson brackets with commu- 
tators multiplied by (7h)-'. We shall consider a Schrédinger representation 
with state function V(p;,..., px) . 

The volume element in the original system is 

do = dx,...d%% 


= 6(Xa941) <e8 6(Xsy4n)dx 1 ane ‘e dx x, 


where 6(x) is the Dirac 6-function. The transformed volume element is 
(29) dv = 8(oy41).-- 5(ou+n)gidpr ...dpx. 


g*(p1, ..., px) is the determinant | dx./dp3 | . Since we want to separate matrix 
elements into two factors depending only on the internal and collective co- 
ordinates respectively, we have imposed the condition on the transformation 
Xa — pg that g* may be written as 


(30) g'(p1, «++, Px) = gu? (p1, wees pw)en' (Past, oc) Pu+n) : 

The operators pi,..., x are Hermitian since they are real functions of 
X1,.-+,Xay and Xy41 =... = Xy4ny = 0. Hermitian representations for the 
momentum operators o1,..., 04, ™w4+1,---, 7x Must contain the derivatives, 
0/d0pq, in the combination 
(31) Da = —1thd/dpa—th /2.dlog {5(a+1) tee 5(bu+n)g*} /Apa. 


The extra term arises since the representations are Hermitian only if the 
relation 


(32) Svo*(DaWa)5(b4041) ee 5(ou+n)gidpr ..-dpx = f (Dae) *Y.5(oau41)-.-dpx 


results from an integration by parts (cf. Pauli 1933, De Witt 1952). The 6- 
functions and their derivatives do not appear in combination with 0/dp, or 
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with those linear forms of 0/dp;, namely 0/0p;—a;;0/0p;, which satisfy 
equation (18) : 
(18) (0/dp:—ai;0/dp;)xn = 0. 


The Poisson bracket expressions (22) and (23) are consistent with represen- 
tations of x, and o; which contain the derivatives in the forms —ihd/dp, and 


—th[d/dp;—(0Xm/Ap:)(Op;/AXm)09/Ap;] , 


respectively. Therefore, the Hermitian representations for these momenta are 


(33) tn = —thd/dpn—th/2.0gn*/dpn 
and 

a; = —th{d/dp;— (Axn/Op;)(0p;/AXn)09/Ap 5] 
(34) +1h/2 . A{ (xn/Op:)(Op;/Axn)o} /Op; 


—ih/2 . [8;;—(A%n/Op:)(8p;/Axn)0]d log gas*/dp;. 


A direct calculation shows that the representation (34) satisfies the quantum 
mechanical form of equation (25), 


(25’) (0p;/AXn)ooj;+0;(Op;/dxXn)o = 0. 


Collective motions are of particular interest for those states in which the 
state function W has the form 


(35) V(p1, sey px) = (puis, sey Pu+n)Xx(P1, sey Pu) ’ 


and for transitions between two such states with the same internal configura- 
tion. Factorizations like (35) are possible in the present formalism since the 
dynamical operators divide into two sets, {p;, o;} and {pm, mn}, such that the 
members of one set commute with each of the members of the other set. The 
state function specifying the internal configuration, x(p1,...,p), is de- 
termined from the state function in the laboratory coordinate representation, 
W[x1,...,X], through the equation 


(36) V[p1, ceey pu] = x(p1, coey pau) o(0, Ho Sg 0) ° 
Hence, knowledge of some of the properties of [x1, . . . , Xx] provides informa- 
tion about x(p1,..., a) . This can be used in calculating the coefficients of 


functions of p, and their derivatives in the dynamical equations of the col- 
lective motions (39). Equation (35) may be rewritten as 


(35’) v= ¢' (putt) Vee y pu+n) ¥[p1, eoey pu] ’ 
where 
(37) ¢’ (pxr+1, sees Pu+n) = (pm+1, corey pu+n){ (0, +a Sy 0)}-". 


Let us consider a state described by a separable state function W and 
characterized by the steady state Schrédinger equation 


(38) H(x1,..., Pu)¥ = EW. 


The Schrédinger equation for the collective motions is 
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(39) Heoulemss, ++» tuen)? = Econ, 
where 
(40) Acou = Sx*Ax8(bar41) vee 5(¢x)gubdpr ...dpy 


for A equal to H and to E. 

The formalism may be extended to take into account the interaction be- 
tween the collective motions of a system whose particles have laboratory 
coordinates x; and other particles with laboratory coordinates y, (a=1,..., Y). 
The coordinates y, and their conjugate momenta gq do not enter into the 
transformation investigated above. The previous discussion applies to the 
extended system by replacing L and H by 

L= Eli aca 5 ties eee 


and 
H = H(p1,..., 4x; Yu-++sQy)- 


CENTER OF MASS AND RELATIVE COORDINATES IN QUANTUM 
MECHANICS 
The simplest example of the formal method presented above is provided by 
the quantum mechanical description of a system of particles in terms of their 
center of mass and relative coordinates. The procedure is simplified by letting 
Xa and pg denote vectors. The operators 0/0xa, 0/dps denote the divergence 
operators V with respect to xq and pg respectively. 


The laboratory coordinates of the M particles, x:,..., x, are related to 
the relative coordinates pi, ..., py and the coordinates of the center of mass, 
pu+i, through the equations 
(41) Xi = pitpmii. 


The internal coordinates satisfy the constraint 

(42) Mipi = 0, 

which prevents the system from undergoing center of mass motion in the 
internal coordinate system. A one-one transformation between the 3M- 
dimensional subspace (42) of the (3A +3)-dimensional p-space and the labora- 
tory configuration space is obtained by introducing a new (identically zero) 
laboratory coordinate x41, 

(43) Xu+1 = mip; = 0. 

The momenta belonging to the coordinates pa, 

(44) Ta = OL(pitpus1, Pit Psi) /Opa, 


satisfy the relation 
(45) Zz T1— Ty+1 = 0, 
t 
which follows from the functional form of the Lagrangian. A new laboratory 
momentum py+41 is defined to be zero and to satisfy the equation 


(46) te = (xp/Opa) pp. 
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The matrix || dxg/dpq || is derived from (41) and (43): 


i *#...# 

ce t..07 1| 

(47) | dx3 / Opa | a 1 . . ° ‘ | 
| ® e . . 

}o 0... 1 Of 

|| 

| mil mol... myl 0| 

' i 


where 1 and 0 are unit and zero 3 X3 matrices. The inverse matrix is calculated 
by the methods of Appendix A: 


(48) |{Ap./dx,!| = 


(1—m,/>-m;,)1 (—m2/dim,)1 ...  (—my/Sm,1 (1/y-m,)1 | 
(—m,/dm,)1 (1—m2/d-m;)1 be (—my/Xm;)1 (1/Xom,)1 
(—ms/Xm)L (= m2/CmyL... (l—mu/Xm) (1/1 

(m,/d-m,)1 (m2/ m,)1 tae (my/dm,)1 (— 1/Ym,)1 


The 7, and pg are therefore related by the equations 
wi = Pitm; Pus, Tus = LPy 
Di = mi—m,/Lm,(Um—P), Puri = (1/Lom,)(L4,—P) = 0. 


(49) 


The Poisson bracket relations between the p,z, 73 variables are 
[pa, pa] = 9, [oi, wi] = (6i;—m;/Lem,)I, 
(50) [pi, Tari] = 9, [oarsi, 75) = (m;/Qom;,)1, 
lexti, Tu+1] = 1, [a, m3] = 0. 


The z; are not independent of the collective coordinate py41 because [ps¢41, 7;] 
does not vanish. However, the linear combinations 


(51) += tee 
have zero Poisson brackets with the collective variables: 
(52) [esrti, 0] = [targ1, oj] = 9. 


We also have 


(53) [p:, 0] = [8i;—m;/Lem,ll, [o:, o,] = 0. 
The momenta o; satisfy the identities 
(54) oe; = 0, 


and are related to the laboratory momenta through the equations 


(55) pi = os tm, /Umyry41. 
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If the system is characterized by the Lagrangian 


(56) L= ¥ 4m, 27+ D v(x1-x,), 


i i<j 
then the Hamiltonian is 


(57 H = Zz oa, /2mi+ wus. /2Dm;+ ¥ v(pi—6))- 


i<j 
For this system, 
(58) o; = Mp; (no summation). 
In the quantized system, pa, oj, and my41 are operators satisfying quantum 
Poisson brackets formally identical to the classical expressions (50), (52), and 
(53). In a Schrédinger representation with state function V(pi, ..., psr41), the 
operators p, are represented by their eigenvalues, and m4: and a; by 


(59) 1 =" 1hd/Opx+1, 
(60) a; = —th{d/dp;—(m;/>m,)>0/dp,). 
The volume element is 
(61) dv = 8(m.p;)(Xm;)'der. . . dows 
= 5(m.p:/2m;) dp... dpy4i. 
For a state with a separable state function 
YW = $(em+1)x(1,-- +» Pa) 
and characterized by the Hamiltonian operator (57), the center of mass motion 
is specified by the Schrédinger equation 


(62) — (hk? /2X0m,) 0°/dpusr’ & = Eco. 


Here 


ae 


(63) E. = J v| et Le 2m. 8/ap.— (m./Xm;) > a/ap.t 
i k 


-> o(o-o1) | x 6(mipi/Xim;) dp... dpm, 


i<j 


and E is the total energy of the system. 
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APPENDIX A. INVERSION OF THE MATRIX ||dx,/dpg| | 
In this appendix only, the transformations are restricted to those having 


the form 
Xi = Xi(Pi, Pmtiy- ++) Px) - 
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With this restriction, the matrix is 











ai; Cinl| 
in Ol 
where 
a; = 8x,/dpi, 
bing = 9%m/ Oj, 
and Cin = 0%;/Opn. 


There is no summation over repeated indices in this appendix, except where it 
is explicitly stated otherwise. This matrix may be factored into the product 


avy, 0 | O43 


643 Cus. 
iO tem hi Mgt 


° | x 0 Dmx 





























where 


/ jo 
Cin = Q; Cin 


and 


, -1 
i — > bmi Cin = — > bn: a Cin- 
t i 


The inverse sania is 








| | 
| 5 iy ae | |x| | - 9 ly - 64; 0 | 
i 0 | — welt 0 bmn 
fai Cim Ann 5,)0;" —a;" :. Cim Amn 
=. mn m 
a Pe Anm Omi a; Line | 
Ld ! 
where 
1! Amn || is the matrix inverse to || Dn || 
and 
oo =, Cm Ann = -—a;'> Cim Amn: 
™m m 


APPENDIX B. EVALUATION OF THE POISSON BRACKETS OF p, AND zg 


The coordinates pz may be expressed as functions of the coordinates 
X1,...,Xx only so that 


[pas Pp] = 0. 
The bracket expression [p,, 7s] may be evaluated by using equation (12) to give 
[Par A] = (Apa/dx;)(Ax;/Ops)[xi, P,] 
= (Opa/dx;)(0x;/Apz) . 


Similarly, 
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(a, ™p] = [(0%;/Opa) pi, (x ;/Aps)P;] 
= (0x;/Apg)(9{ Ox ;/Opa} /x;)pi—(Ixi/Apa)(9{ Ax ;/App} /Ax i) p; 
= (dx;/Opg)(0{dx;/Apa} /Ax;)(Apy/Ax ;) ry 
— (0x ;/Opa)(9{ 0x ;/Apg} /Axi)(Apy/Ix;)ry. 


The bracket expressions (22) and (23) are evaluated as follows: 
[pay ®m) = (Opa/dxi)(Ox;/Apm) = (Apa /Axg) (Axp/Ipm) 


Sam; 

(Pn, o1] = (Opn/Ox,) (0x4 /Op;) — (Opn/AX%) (Ox 4%/Ap;) (89; /8Xm)0(OXm/Op:) 

— (Opn /OXm) (OXm/Op:) + (Apn/AXp) (AXp/Ap;) (Ap s/AXm) 0(AXm/ Api) 
= 0, 

(pi, oj] = (Opi/Ox~)(Oxx/Ap;) — (OXm/Op;) (Apr/AXm)0(Opi/Ix x) (Ox%/Ipx) 
= 63;—(0p;/AXn)(0Xn/Op;) — (0Xm/Ap;)(Opi/AXm)o 

+ (0%Xm/Ap;)(Opn/OXm) 0(Opi/IXn) (OXn/ Opn) 

= 5:;—(0Xm/Op;)(Opi/IXm)o, 

(an, Tm] = (0x ;/Apm)(O{ Ox :/Apn} /Ox;)pi— (Ox i/Opn)(9{ Ix ;/Ipm} /Axi)p; 
= (0°x;/dpmOpn)Pi— (0°x;/ApnOpm) pj 
= 0, 

(dx ;/Ap;)(8{Ax%/Apn} /Ax;)px—(Ax;/Apn)(O{Oxz/dp;} /dx;) pe 

— (0Xm/Ap:)(A{0x%/Apn} /OXm) Pky 

[ny 71) = [—(0Xm/Op:) + (8Xm/Ap;) (Op j/AXp) 0(AXp/Api)](A{AxX%/Apn} /AXm) Px 
= 0. 


A straightforward, but lengthy, calculation using 
[wi, 73] = —(0Xn/Ap;)(A{Ox%/Ap:} /AXn) Pet (AXn/Op;)(A{Ax4%/Ap;} /AXn) Pe 
gives the value of [¢;, oj] shown in equations (28). 
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ENCLOSED QUANTUM MECHANICAL SYSTEMS! 


By T. E. HuLt anp R. S. JuLius 


ABSTRACT 
A brief description is given of the eigenvalue problems associated with enclosed 
quantum mechanical systems and of some attempts to deal with these problems. 
Another method is developed which leads to a general asymptotic formula for 
the eigenvalues. This formula yields a simple asymptotic approximation to 
the eigenvalue in each particular case, once the eigenfunction of the corre- 
sponding unrestricted system is known. 
INTRODUCTION 
The type of problem to be considered can be described by first considering 
a simple example. The differential equation for the linear oscillator is 
9 
du 2 
—3z—xutru = 0. 
dx 
Usually we want to find the solutions of this equation which vanish at the 
singular points x = +o. We then obtain the eigenvalues X, and the corre- 
sponding normalized eigenfunctions u, where 


ho = 1, Ug = wi tlte-2? 22, 
A = 3, (4/x)"/*xe-#*2, 


Ae = 3; U2 (4r)—!/4(2x? — 1)e-7? 2, 


2 
ll 


However, for some applications we might consider the oscillator to be 
confined to an enclosure with walls at x = +a. (For example, we might 
approximate to the effect of pressure in this way.) We then have to find 
solutions of the above equation which vanish at x = ka. 

For certain values of a, besides a = ©, we already know an exact solution. 
For example, since “2 vanishes at x = +1/+/2, then uw represents the ground 
state of the oscillator when a = 1/+/2—except that ue is no longer properly 
normalized. When a = 1/+/2 the eigenvalue for the ground state is therefore 5. 
Determining other such results enables us to draw a graph of \ against a for 
each state (Fig. 1). The eigenvalues for any particular value of a can now be 
determined approximately from these graphs. 

Similar graphs can be drawn for the three-dimensional oscillator or the 
hydrogen atom confined to a sphere, or for the restricted rigid rotator. Their 
graphs are more complicated because their excited states depend on two 
quantum numbers instead of one. 

Better approximations to the eigenvalues can be obtained, in general, only 
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First Excited State 





Ground State 


' 2 3 + 


be a 


Fic. 1. Dependence of eigenvalue \ on size of enclosure containing a linear oscillator. 


by numerical methods. But analytical expressions for \ can be obtained in 
special situations. Such expressions are usually asymptotic representations of X. 

The earliest result of this kind appears in a paper by Michels, de Boer, and 
Bijl (1937), who consider the ground state of the enclosed hydrogen atom. 
Their approximation improves as the radius of the bounding sphere increases, 
but it then becomes less manageable. De Groot and ten Seldam (1946) remove 
this difficulty and apply the method to 2s and 2p states and also give numerical 
results for radii which are not large. Analytical results for all states and dif- 
ferent sizes of radii are given by Dingle (1953). Related results are obtained 
by ter Haar (1946) for another quantum mechanical problem, and by 
Abramowitz (1953) for a heat transfer problem. 

Chandrasekhar (1943) uses the method of Michels ef al. to deal with the 
first excited state of the enclosed linear oscillator in connection with an astro- 
physical problem. Auluck and Kothari (1945) find an asymptotic expression for 
the energy of the oscillator which is valid for small a; their expression for large 
a is not valid because of an incorrect use of asymptotic expansions. 

Sommerfeld and Hartmann (1940) consider the restricted rotator. They 
give the graphs corresponding to Fig. 1 and they analyze in detail the limiting 
case when the rotator is only slightly restricted. Related results are given by 
Siegel et al. (1952) in connection with a flow problem. 

Further discussion of the above problems is given by Cotrell (1951), and 
ten Seldam and de Groot (1952). 

We proceed now to develop a method which leads to a general asymptotic 
formula for the eigenvalues of an enclosed quantum mechanical system. The 
formula is valid in the ‘‘neighborhood”’ of any known solution. With the linear 
oscillator, for example, the method is valid near a = & or, in the ground state, 
near a = 1/,/2. In the last section we give simple asymptotic approximations 








916 CANADIAN JOURNAL OF PHYSICS. VOL. 34, 1956 


to the formula for each of the problems mentioned above. These approxi- 
mations are given only for the cases in which the known solution is the solution 
of the corresponding unrestricted problem. 


GENERAL FORMULA 


We suppose that we already know an eigenvalue \ and a corresponding 
normalized eigenfunction u for the problem whose differential equation is 
(1) u'+r(x)ut+trAu = 0 
and whose boundary conditions are u(0) = 0 (or u’(0) = 0) and u(d) = 0. 
Here r(x) is assumed to be continuous for0 < x < } except that the origin may 
be a regular singular point of the differential equation. If d is finite, b is assumed 
to be an ordinary point or a regular singular point of the differential equation. 

We now have to consider the “perturbed” problem whose differential 
equation is 

y" +r(x)y+(A+Ad)y = 0 

and whose boundary conditions are the same as those for the‘‘unperturbed”’ 
problem, except that 5 is replaced by a. Here r(x) must be continuous for 
0 <x <asothata < dif d isa singular point. 

Suppose that v is any solution of (1) such that the Wronskian wv’ —u’v = 1, 
and-‘suppose we normalize y by requiring that y(0) = u(0), y’(0) = u’(0). 
Then it is easily verified that 


y(x) = u(x) —Ar J (0x) u(t) u(x) 9) 9) 


Putting x = a in this equation and requiring that y(a) = 0 we obtain 


(2) A\ = u(a) / (oa) Ju y(t) dt —u(a) Jo y(t) it) : 


To find what happens to this expression for AX as a > 6 we first note that 
y(x) depends continuously on Ad; therefore y(x) — u(x) as AA — 0, and hence 
as a — b. If we also assume for the moment that 


(3) ly~)|<M, O<x <a, 
where M is independent of a if a is sufficiently close to b, then we have 


(4) v(a) Ju y(t) dt~v(a), ab, 


where we have used the fact that u(x) is normalized. 
To deal with the second term in the denominator of (2) we first write 


(5) |ue) foo uy < |ula) J oH 9 at] +} u(a) fo 9 a 














< a)|(atar fiat), 


where A is some finite constant. If c is chosen so that u(t) has no zero in the 
interval c << t < b (orb <t <cifa> b), wecan, in this interval, put 
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(6) v(t) = u(t) fro dt. 


The behavior as a — 6 of the right side of (5) is therefore governed by the 
behavior of u(t) as t > 6. Only two possibilities need be considered : 


(i) Ifb << @, 
u(t) ~~ B(t—d)2, t— b, 
where B, a are constants and a > 1/2. From (6) we then obtain 
—p)! yee = 
ee (tb)! In(t—b)/B, a = 4,48, 
(t—b)**/B(1—2a), a>toob. 
(ii) fb =o, 
u(t) ~ Btze8*, t—b, 
where B, a, 8, y are constants and 8, y > 0. From (6) we then obtain 
v(t) ~ &"/2BByit, tb. 
Substituting from either of these possibilities into (5) leads to an asymptotic 
expression which is negligible in comparison with the expression in (4). The 
second term in the denominator of (2) can therefore be neglected as a — bd. 


Using this result, along with (4) and (6), we finally obtain from (2) the asymp- 
totic formula 


(7) A\ ~ i/ fiero dt, a—b. 


We have still to establish (3). This can be done for each case separately; we 
shall consider one typical case: the even states of the linear oscillator with 
b =o. Here 

yl!’ —x?y+(2n+1+Ar)y = 0 


for some even integer m. For the solutions we obtain 
@o 
—z2/2 2k 
y(x) = "DE bx 
k=0 


so that y’(0) = u’(0) = 0, and y(0) = u(0) if bo is properly chosen, and where 
bear = [(42 —2n—Ad)/(2k+2)(2k+1)]dy. If a is sufficiently large the AA 
which makes y(a) = 0 will certainly satisfy | AN | < 4. Then 














= bax" < (x /a)"*? yi ba , O<xK<a, 
n/2+1 m/2+1 
n/2 
= (x/a)"**?| 3° ba™| , since y(a) = 0, 
0 








< (x/ay**¥ [onl a” 


n/2 
< dX |de| 2 


n/2 


< = by att 
0 
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where | by | < 5,*, and 5,* is independent of Ad and hence of a. We therefore 
have 


n/2 
ly(x)| < 26°"? > b,*x" < M 
0 


for 0 < x < aand where M is independent of a, as required for (3). 


SPECIAL RESULTS 


Except in special cases the general formula (7) is not easily evaluated. 
However, an asymptotic approximation to (7) is easily obtained from the 
expression for u(t). If 6 is a singular point these expressions are well known 
and we shall list the corresponding results for each of the problems mentioned 


earlier. 
The method of calculation is illustrated with the linear oscillator problem 
where (Morse and Feshbach (1953)) 


\ —n/2 2 mE paca 2 <a — 12 
u(t) ee ate (1) eG (e*). 


The normalization has been adjusted so that fj r(t)dt = 1, as required for (7). 


We obtain 
u(t) Sai Qri2+ii2( 9 VN —V2g—l/4pne— tie too, 


so that 


A\~ armies / | pas ae 


+2 \—1 .—1/2 n+] —a? 
we al) ee”: 8: fl ee. 


For the linear oscillator we therefore have 
y"—xy+(A+Ar\)y = 0, y(—a) = x(a) = 0, 
A = 2n+1, TO) ee he cn: 
AN ~ 2%42(n !)—le 1 2gintig—2 ao, 
For the three-dimensional oscillator we obtain, for / = 0,1, 2,..., 
y” —(P+114+1)/x*)y+(A+ArA)y = 0, y(0) = y(a) = 0, 
A = 4n—2/+4+3, n= 1,1+-1,1+2,..., 
AN ~ 4a4n-2!+3¢e—9° /(m —])! P(n+3/2), ao, 
For the hydrogen atom we obtain, for] = 0, 1, 2,..., 
yy” +(2/x—l(1+1)/x*)y+(A+AA)y = 0, (0) = y(a) = 0,7 
A= —1/n’, n=1+1,/+2,/4+3,..., 
Ad ~ 22% +1g2ne-24/n /42n+3(y —]—1)! (n+1)!, ao, 
For the rigid rotator we obtain, for m = 0, 1, 2,...and for e small and positive, 
y’’ —((m? —1/4)/sin’x+1/4)y+(A+4A)y = 0, =e) = v(x) = 0, 
A = LU(l+1), 1 =m, m+1, m+2,.... 
If m > 0, AX~ m(1+m)!(214+1)2"/2?"—(1—m) !(m !)?, e— 0. 


If m = 0, AN~ —(2/+1)/Ine, e— 0. 
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(The expression for m = 0 is twice the one obtained by Sommerfeld and 
Hartmann (1940). This is to be expected since, in our notation, they used the 
boundary conditions y(0) = y(r—e) = 0.) 


To describe the accuracy of the above results let us consider the extreme 
case of a hydrogen atom in a sphere of radius five, i.e., five times the Bohr 
radius, which corresponds to a pressure of nearly 5000 atmospheres. In this 
case formula (7) can be evaluated in terms of the exponential integral. The 
result is equivalent to that given by de Groot and ten Seldam (1946) and gives 
AX to within 1% and A+AdA to within 0.01%. On the other hand, the simple 
approximation found in this section gives AA to within 30% and A+AA to 
within 0.3%. 
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AN OPTIMUM PARAMETER BETA-RAY SPECTROMETER! 


By C. GEOFFRION AND G. GIROUX? 


ABSTRACT 


Un spectrométre beta semi-circulaire a été construit en donnant aux divers 
parameétres de l'appareil des valeurs telles que l’intensité maximum d’une raie 
d’électrons monoénergétiques soit la plus grande possible 4 une résolution 
donnée. L’appareil balaie un spectre de fagon automatique, mesure des raies 
avec une limite de résolution atteignant 0.07% et peut détecter des électrons 
d’énergie variant entre 0 et 4 Mev. 


INTRODUCTION 


During the last decade, many studies relative to the improvement of beta- 
ray spectrometers have been published. One of us (Geoffrion 1949) has in 
particular studied the optimum conditions of operation of the 180° instrument. 
A spectrometer has been built along the lines suggested by this study. The 
apparatus has high resolution, is automatically operated, and is equipped 
with an electrostatic accelerator for the detection of very low energy electrons. 
It is the purpose of this paper to describe the spectrometer and its principal 
features. 


OPTIMUM CONDITIONS OF OPERATION 


According to the above-mentioned paper, there exists for the 180° beta-ray 
spectrometer a particular choice of the various parameters of the apparatus, 
such as source and slits dimensions, that will provide a beam of monoenergetic 
electrons of maximum intensity for a given limit of resolution. These con- 
ditions also provide maximum intensity beams from continuous spectra for a 
given limit of resolution. These conditions are: (1) the width of the source 
and of the exit slit should both be equal to (4r/3)(Ap/p); (2) the 90°-position 
diaphragm width should be 2r(2/3)*(Ap/p)?; (3) the height of the source, 
the diaphragm, and the exit slit should be 227(1/3)?(Ap/p)!. In these expres- 
sions, r is the mean radius of curvature of the electron trajectories and Ap/p 
is the limit of resolution, which throughout this paper refers to the limit 
measured at half the peak height. 

Under these conditions, the relation between the maximum intensity of a 
line and the limit of resolution is: 


te = 0.457 Kr? (Ap/p)5”, 


where K is the number of monoenergetic electrons emitted by the source in 
all directions per unit area, per unit time. In the case of a continuous spectrum, 
the spectral density, i.e. the number of electrons of momentum between p 
and p+Ap reaching the detector per unit time divided by Ap, is then given 
by 
W = 0.207 kr?(Ap/p)5”, 
1Manuscript received April 23, 1956. 
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where & is the number of electrons of a continuous spectrum emitted by the 
source in all directions per unit area, per unit time and per unit momentum 
interval. 

THE SPECTROMETER CHAMBER 


In our optimum-parameter 180° spectrometer, the radius of curvature of 
the electron trajectories is 30.5 cm. and the different parameters have been 
chosen so as to provide, for normal operation, a limit of resolution of 0.25%. 
According to the previously-mentioned optimum conditions, the different 
parameters of the spectrometer should then be: (1) width of the source and 
exit slit, 1 mm.; (2) width of the diaphragm, 2.4 cm.; (3) height of the source, 
diaphragm, and exit slit, 5.34 cm. 

The spectrometer chamber is shown in Fig. 1. It is made of brass and is lead 
soldered; the three ports give access to source, diaphragm, and exit slit, and 
have vacuum-tight metal-to-metal (with vacuum grease) or O-ring joints. 
The chamber is evacuated with a 50 liter/min. diffusion pump backed by a 
mechanical one, and a vacuum of the order of 10-5 mm. Hg can easily be 
obtained. The source holder, the diaphragm, and the exit slit are made of 
1/4 in. thick aluminum plates, sliding on grooves to reproducible positions. 
Thin auxiliary Al fins are disposed along the trajectory limits so as to baffle 
and reduce scattered electrons. The Geiger counter is rigidly attached to the 
exit slit door. 


THE MAGNETIC FIELD 


The homogeneous field necessary for the spectrometer is produced by an 
electromagnet of the double-mushroom or yo-yo type (Figs. 1 and 2). It is 





— 12 in— 


Fic. 1. Spectrometer magnet and chamber. O, magnet core. C, magnet coil. d, chamber 
ports. S, source holder. D, 90°-position diaphragm. E, exit slit. G, detector holder. V, to 
vacuum system. B, electron beam. 

Fic. 2. Magnet cross-section. A, air gap. O, core. C. coil. 
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made of ordinary cast iron machined to a tolerance of 0.001 in. and its field is 
sufficiently homogeneous without shimming in the region of the air-gap 
utilized; any fringing effects present will improve the focusing. The coil 
consists of 30,000 turns of Formex 24 wire; it is fed by a stabilized current 
power supply and, at maximum field (1200 gauss), dissipates 24 watts. However, 
the magnet is always used below saturation of any section, and 3 Mev. electrons 
are focused by a 400 gauss field produced by 100 ma. 

The strength of the magnet field is determined and stabilized by comparison 
with the magnetic field of a pair of Helmholtz coils. Two identical search coils 
rotating at 30 cycles/sec. (along one of their diameters) in each of the spectrom- 
eter and Helmholtz fields generate electromotive forces which are compared 
in a given ratio through an attenuator. A null difference between the two 
attenuated signals indicates that the two fields are in the chosen ratio. The 
difference, if any, after power amplification served to actuate a servomechanism 
which controls the current flowing in the magnet. The magnet field is thus 
determined and stabilized by the choice of the attenuation factor and by the 
strength of the Helmholtz field. 

The two rotating coils have a mean diameter of 1 cm. and contain 4300 turns 
of No. 34 copper wire; they are fixed at the ends of a stainless steel rod, 1 in. 
in diameter and 3 ft. long. The rod is horizontally supported by ball bearings 
and rotates at 30 cycles/sec. along its axis, driven by a synchronous motor. 
Three rotating contacts transmit to the attenuator the electromotive forces 
generated by the coils. Each contact system (Fig. 3) is made of a brass pulley 
attached to the steel rod but insulated from it and connected to one coil 





Fic. 3. Rotating contact system. 
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terminal, and of a second brass pulley fixed to a brass shaft, free to rotate on 
ball bearings in a vertical position. At the extremity of the shaft is a platinum 
needle dipping in a cup of mercury. The needle shaft is driven by a steel 
spring belt joining the two pulleys and the signal is picked up at the mercury 
cup with a platinum wire. In order to avoid mercury oxidation, a small amount 
of silicone oil covers the mercury surface. This system has a few hundredths 
of an ohm of resistance which has remained constant, even after hundreds of 
hours of operation. It is far better than any other contact system tried. 

Through the rotating contacts, one terminal of each coil is grounded and the 
others are connected to the attenuator (total resistance 100 kQ in 10 Q steps); 
the coils are arranged in such a way that their signals are out of phase. From 
the attenuator, the difference or ‘‘error signal’’ passes to an amplifier, con- 
sisting of two narrow band stages (Sturvenant 1947) plus an additional 
amplifier stage; the over-all gain is around 10°. After power amplification, 
the signal feeds the secondary of a two-phase reversible motor; the motor 
primary is fed by a power amplified signal coming from a photocell which 
receives a light beam chopped by reflecting and non-reflecting parts of the 
rotating steel rod. The primary voltage is thus stable in amplitude (100 v.), 
phase, and frequency; the motor will turn clockwise or counterclockwise 
according to the relative phase (+90°) of the two voltages. This phase, in 
turn, is determined by a greater or a smaller value of the magnet field relative 
to the Helmholtz coil field. The servomechanism control is completed by 
having the motor drive a 40-turn Helipot potentiometer, the sliding contact 
of which fixes the grid polarity of two 6AG7 tubes controlling the magnet 
current. The power supply also contains a conventional voltage stabilizer. 

In order to avoid oscillation of the system, due to the magnetic inertia of 
the magnet, a d-c. generator was coupled to the motor and an appropriate 
fraction of the generated voltage is added to or subtracted from the grid 
polarity of the 6AG7’s, making an overcorrection that disappears when null 
difference is reached. The error signal amplifier is also equipped with an 
automatic gain control so that large error signals do not saturate it. Figure 
4 illustrates the stabilizing system. 





Fic. 4. Field stabilizing system. M, magnet. H, Helmholtz coils. R, rotating steel rod. 
c, rotating search coils. v, rotating contacts. n, light chopping mirror. }, light bulb. p, photocell. 
a, attenuator. Aj, error signal amplifier. As, stable signal amplifier. m1, 30 c.p.s. two-phase 
reversible motor. g, d-c. generator. S,, magnet power supply. u, converter. F, reference voltage. 
P, potentiometer. A;, error signal amplifier. m2, 60 c.p.s. two-phase reversible motor. S2, 
Helmholtz coils power supply. 
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The Helmholtz coils (Fig. 5) have a mean radius of 11.7 cm. and each 
contains 8000 turns of No. 27 copper wire; each coil has a conical form, so as 
to satisfy the Helmholtz condition for each turn of wire and ensure a field 
uniform to at least a few parts in ten thousands in the central region where 
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FIG. §. Helmholtz coils. H, large Helmholtz coils. C, coils. c, cradle. T, small Helmholtz 
coils. @ and 8, pivots. 


the search coil rotates. The field generated with 100 ma. in the coils is 60 gauss. 
The Helmholtz coils can be rotated about the axis passing through the two 
rotating coils so as to adjust the relative phase between the two generated 
signals. In order to avoid any error due to the presence of the earth’s magnetic 
field, a second pair of smaller Helmholtz coils, concentric to the first one, 
cancels this field in the search coil region. 

The main Helmholtz field current is stabilized and fixed by servomechanism 
techniques by comparison with a stable reference voltage (Fig. 4). The current 
flowing through the coils also produces a voltage drop in a standard manganin 
resistor; this voltage drop is compared to the reference voltage by a 60 cycle 
L. and N. converter. The resulting error signal serves, after amplification, to 
actuate a two-phase reversible motor the primary of which is actuated by the 
110 v. line. The motor drives a 40-turn Helipot, the sliding contact of which 
fixes the grid polarity of 6AG7 valves controlling the coils current. The rest 
of the power supply is a conventional voltage stabilizer. 

The reference voltage is a fraction of an 11-volt stack of mercury batteries. 
The voltage divider (Fig. 6) is made of a 10 kQ resistor variable in steps of 
one ohm, and of two 50 kQ continuously variable Helipot potentiometers. 
The two Helipots serve to determine the general operating region and the 
precise voltage reference value is given by the four controls of the step 
divider. 

A second standard manganin resistor (10 ohms), in series with the Helmholtz 
coils, permits accurate potentiometer measurements of the coil current; 
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Fic. 6. Reference voltage divider. V, 11 v. mercury battery. H, 50 K variable Helipot 
resistors. Ri, 10 4 10-step variable resistor. Re, 1 K, 10-step variable resistor. R3, 100 ohm, 
10-step variable resistors. R4, 10 ohm, 10-step variable resistor. 


J 





if the whole system has been calibrated, these current measurements give the 
Hr value of the registered electrons. 
THE DETECTOR 

The detector used is a Geiger—-Muller counter, placed with its axis parallel 
to the exit slit of the spectrometer. It has a very thin Zapon window and is 
equipped with an electrostatic accelerator allowing detection of very slow 
electrons down to zero energy with a few transmission corrections. This 
accelerator—detector has been described previously (Giroux and Geoffrion 
1956). 

The pulses from the counter are sent to a preamplifier and then to a scaler. 
The scaler includes three decimal counting units and one preset decimal 
counting unit; with a suitable combination of the preset unit and of one, two, 
or three other units, the system becomes a scaler of any number from 1 to 
10, multiplied by 10 to any power from 0 to 3. The output pulse from the 
scaler serves to actuate a time-recorder (Streeter-Amet) which counts and 
prints the time necessary for the counter to detect the number of electrons 
selected by the scaler setting. 

AUTOMATIC OPERATION 

Automatic operation of the spectrometer has been achieved in the following 
way: first, the scaler and the time-recorder have been modified so that they 
reset themselves to zero at the end of each counting period. Second, the output 
pulse from the scaler, which actuates the time-recorder, is also used to change 
the magnet field by a preset increment and to start a delay-timer. At the end 
of the delay period (time necessary for the magnetic field to assume a new 
value), the scaler is automatically put in operation. This cycle can be repeated 
again and again, and a whole spectrum is then scanned. 

Since the magnet field value depends on the Helmholtz field which in its 
turn is fixed by the reference voltage already mentioned, the automatic 
operation needs only to control the voltage divider giving this reference 
voltage. The divider consists of two continuously variable potentiometers, 
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fixing the range of operation, and a 10 kQ step divider with four controls 
determining the units, the tens, the hundreds, and the thousands. Using 
relay technique, each control was arranged so that an incoming pulse increases 
its setting by one and so that when it reaches 10, it resets to zero and sends a 
pulse to the next control relay. The incoming pulses from the scaler may be 
sent to any of the control relays, and, for versatility, a mechanism permits 
the multiplication of an incoming pulse by any preset number from 1 to 10. 
The magnet field can thus be increased by any preset value. The system may 
also be arranged so that the field is decreased by the preset amount. In order 
to know the Hr values corresponding to reference voltages, it is necessary 
to make for some of these voltages a potentiometric measurement of the 








Fic. 7. Automatic operation system. M, magnet. R, rotating rod. H, Helmholtz coils. 


G, counter. T, high tension source. S, scaler. Sc, scaler control. n, ‘‘ott’’ signal. 0, ‘‘on’’ signal. 


C, printer-chronometer. p, “off” signal. s, ‘‘on”’ signal. D, delay-timer. 7, ‘‘on”’ signal. u, ‘‘off”’ 
signal. F, reference voltage relay system. Hc, Helmholtz coils current control. 


Helmholtz coils’ current. A graph of Hr vs. reference voltage serves to translate 
all reference voltage values into Hr values. Fig. 7 shows the automatic operating 
system. 

BEHAVIOR AND RESULTS 

To verify the precision, stability, and resolving power of the spectrometer, 
a series of trial measurements has been made. 

Calibration of the apparatus has been made using the F and A lines of 
thorium active deposit. Using the Hr values given to 1 part in 10,000 by 
Siegbahn (1955), it was possible, under optimum conditions, to measure the 
Hr values of the B, I, and J lines agreeing within 1 part in 3000 with the 
Siegbahn values. Under normal conditions of operation, the apparatus is 
able to determine Hr line values within 1 part in 1000. 

With the F line, the limit of resolution of 0.25%, predicted for the slit and 
source dimensions previously mentioned, has been verified as shown on Fig. 8. 
In this case the successive increments of the magnetic field were 1 part in 
2000. This is the usual mode of operation of the spectrometer. 

Figures 9 and 10 represent the F and I lines with a higher resolution. Slits 
and source dimensions have been modified so as to give a limit of resolution 
of 0.06%. The two lines have widths of 0.103% and 0.07,% respectively; 
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Fic. 8. Thorium F-line at theoretical 0.25% limit of resolution. 

Fic. 9. Thorium F-line at theoretical 0.06% limit of resolution. 

Fic. 10. Thorium I-line at theoretical 0.06% limit of resolution. 

Fic. 11. Thermo-electrons accelerated to 152 ev. The resolution is poor owing to initial 
electron thermal energy spread. 
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the fact that the F-line is larger can be explained by taking into account the 
natural line width as already observed by Siegbahn (1955); the increase in 
line width is here the same as that observed by this author. For these measure- 
ments, the relative automatic magnetic field increments have been reduced to 
1 part in 13,000. 

To verify the ability of the spectrometer to register low energy electrons, 
an incandescent filament has been used as source and the thermal electrons 
emitted have been accelerated at the source by 152 ev. The measured line is 
shown in Fig. 11. 


CONCLUSION 


As can be seen from the different tests made with the spectrometer, the 
apparatus has good resolution and stability even with automatic operation 
and is able to measure easily electrons in the very low energy range. The 
measured maximum intensity of lines has also verified that its luminosity is in 
perfect agreement with the theory. Unfortunately, the results do not permit a 
comparison between the usual type spectrometer (non-optimum conditions) 
and this optimum parameter one. Theory however predicts a gain of a factor 
5 in favor of the latter with respect to maximum line intensity. 

The transmission of the apparatus is 0.075% at a limit of resolution of 
0.25%, and 0.019% at a limit of resolution of 0.06% which is smaller than in 
double-focusing or lens spectrometers. However this is compensated by a 
larger source area which in this case is 0:5 and 0.07 cm.?, since the number of 
electrons reaching the detector depends both on transmission and on source 
area, 

In conclusion, we believe that this spectrometer has a yield almost equal 
to some existing double-focusing spectrometers of the same radius. However, 
the latter would always be more luminous for the same resolving power and 
comparable magnet size if appropriate optimum conditions are applied. 
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STANDARD DEVIATION OF DEAD TIME CORRECTION 
IN COUNTERS! 


By L. LoRNE CAMPBELL 


ABSTRACT 


Because of the finite recovery time of counters a correction must be applied to 
the number of registrations which the counter makes so that a better approxima- 
tion to the true rate at which events occur may be obtained. It is shown here 
that, for random events, the standard deviation of the corrected rate is approxi- 
mately [v(1-+»7)/t]4, where v is the mean rate at which events occur, 7 is the 
recovery time of the counter, and ¢ is the duration of the observation. This value 
of the standard deviation is an asymptotic value which holds for large values of ¢. 
As an intermediate result, an asymptotic expression is obtained for the value of 
the mth moment of the number of registrations in time ¢. The method used here 
could also be applied to determine asymptotic expressions for the expected 
values of other functions of the number of registrations. 


I. INTRODUCTION 


Because of the resolving time of counters it is often necessary to correct the 
observed number of counts to allow for the number of events which occur 
while the counter is locked owing to the registration of previous events. A 
similar correction is necessary when events are recorded on charts by pen 
recorders because of the width of the line on the chart. This ‘dead time’ cor- 
rection has been treated by Kurbatov and Mann (1945), Feller (1948), 
Ramakrishnan and Mathews (1953), and Hull and Wolfe (1954). The object 
of the present paper is to estimate the standard deviation of the counts which 
are obtained by applying the dead time correction. The estimate is in the form 
of an asymptotic expression, which is useful when the observation time is large. 

It is assumed that events occur randomly and independently at the rate v. 
That is, the probability that an event will occur in a small interval of time, of 
duration At, is approximately vAt, independently of how many events precede 
or succeed this event. Then the number of events which occur in a time ¢ is 
governed by the Poisson distribution. The probability that exactly events 
will occur in an interval of duration ¢ is 


(1) Pnalt = ir er. 


It is also assumed that the resolving time of the counter is non-extended. 
This means that when the counter registers an event it is locked for a fixed 
time, 7, after which it resumes its ability to register events. Thus, the type of 
counter considered here can register an event provided it has not registered an 
event during the interval T preceding this event. This is the Type I counter of 
Feller (1948) and of Ramakrishnan and Mathews (1953). The type of counter 
in which the recovery time may be extended by subsequent events while the 
counter is locked will not be considered here. 


1Manuscript received May 25, 1956. 
Contribution from the Radio Physics Laboratory, Defence Research Board, Ottawa, 
Canada. The work was performed under project PCC No. D48-28-35-05. 
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II. DETERMINATION OF THE MOMENTS OF THE OBSERVED COUNTS 


Let an event be registered at the beginning of an interval of duration ¢ and 
let P,(t) be the probability that exactly subsequent events are registered 
during the interval ¢. Then Ramakrishnan and Mathews show that 


_ —¥(t—T) 
(2) P| =° o> 2), 
said (t<T), 


and 


¢ ? a y(n, v[t—nT]) Wiens! 
(3) Prt) = "Gay UAT) 


ee iP U(t—[n+1]T) (nm = 1,2,...), 


Nn. 


where U(x) is the Heaviside unit function, given by 


(4) v0)" eons 
= @ (x < 0), 


and y(n, x) is the incomplete gamma function, given by 
4 

(5) y(n, x) = J e ‘t” "dt. 
0 


When 1 is an integer, y(, x) is given by 
> y(n, x) agin a" 
YAM AT _ y- x. we, 0: 
(6) (n—1)! om mao m! - 2, ) 


The probabilities given by equations (2) and (3) differ slightly from those 
which would be obtained from the results of Hull and Wolfe because it has 
been assumed that an event is registered at tf = 0. When the dead time, 7, 
vanishes, it is not difficult to show that P,,(¢t) reduces to the Poisson distri- 
bution p,(¢). 

The specific aim of this section is to obtain an asymptotic representation of 
n™ for large values of t, where n™ is the mth moment of the number of regis- 


trations in the time f. That is, 2™ is given by 
(7) n™(t) = Do n™P,(t). 
n=0 


Only positive integer values of m are treated here. Clearly n%(t) = 1. The 
asymptotic representation is easily obtained by the use of the Laplace trans- 
form. The Laplace transform of a function f(t) will be written Af(t), so that 


(8) L f(t) = J e'f(t)dt. 
0 
It is shown by Ramakrishnan and Mathews (1953) that 
n —npT n+1. —(n+l)pT 
9 LP, (t Ss ts ence ae DOs a ie 
ss O-5oby por ' 


Hence 
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a p+ y—pe 97 a PE 
(10) n'(t) = orca » G+ (m = 1,2,...), 
while 
(11) YF nt) = 1/p. 
It is not difficult to show that 
ae —pT m—1 
(12) FE ney: = pr De mCe Ln (t) (m = 1,2,...), 


ce 
where ,,C; is the binomial coefficient. Hence Y’n™ can be expressed in terms of 
Sn* fork = 0,1,...,m—1. 
An asymptotic expression for 7™(t) can now be obtained with the aid of a 
Tauberian theorem given by Doetsch (1943) and van der Pol and Bremmer 
(1950). In its simplest form this theorem states that 


(13) lim h(t) = lim gr". h(t), 


m ! 
ce Mt 5.0 


provided that the limit on the right exists and that A(t) is monotonic for f > 0. 
Since n”(t) is the expected value of the mth power of the number of regis- 
trations in time /¢, it is clearly a monotonic non-decreasing function of f¢. 
Therefore, 

mM 


(14) Sas _ = lim gre iY % n™(t). 


tra m! » 50 


It follows, by induction, from equation (12) that 


ss . nn(t) . ¥ 
15) l #.{ ak 
( ) ~~ t I+v7 


Similarly, it can be shown by induction that 


— y m m—1 9 
(16) im "| wR —( 4) i sine pyaar [m—(1+vT)’]. 


The derivation of equation (16) is complicated by the fact that the function 


a4 ‘*. 7 
ie I+1T 


is not necessarily monotonic for ¢ > 0. It will be demonstrated in the Appendix, 
however, that equation (13) may still be used for this function. An interesting 
result of this demonstration is that a lower bound can be found for n™(t). The 
bound is given by 


i om ie 
(17) n > (28 (> Tf). 

From equations (15) and (16), the asymptotic representation 

W(t) = (—#4_)" soles O-ta m1 g(gn-t 


follows. In equation (18), o(t"-') is a function of smaller order than #”™~'. 
That is, 
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m—1 
(19) lim 2° = 0. 


too 
The cases m = 1 and m = 2 yield 
- vt (2vT+v°T") 


i = Ta” shat 


x2 


and 
= ve \?  (1—2vT—»*T"*) , 
21 Pe Soren = a ae t t). 
ey . (245) Tr 
The quantity # is the expected number of registrations in the time ¢. The 
second term on the right of equation (20) differs from the results given by 
Feller (1948) and Hull and Wolfe (1954) because of the slightly different 
formulation of the problem. If t is replaced by t+T in the right-hand side of 
equation (20) the result is identical with that of Feller and of Hull and Wolfe. 
The standard deviation of the number of registrations, ¢,, is given by 





* = Gantt: 
Thus, as pointed out by Feller, one may not expect % and ¢,? to be equal unless 
T vanishes. 

It is usual to assume that ¢ is so large that only the first term on the right- 
hand side of equation (20) is significant. If this is assumed, and equation (20) 
is solved for vy, the result is 
(23) y = n/(t—AT). 

Equation (23) gives the usual dead time correction. 

The value of x” for smaller values of ¢ may be obtained directly from 
equations (2), (3), and (7). Hull and Wolfe (1954) give an alternative expression 
for % which is more useful for numerical work. 


III. MEAN AND STANDARD DEVIATION OF THE CORRECTED VALUES 


Let No be the rate of registration of events. Thus, 


(24) No = n/t 
in the notation used here. Then the corrected rate, N,, is defined by 
(25) N. = No/(l1—NoT) . 


From equation (23) it is seen that if No = 7/t, then N, = ». 
In this section asymptotic representations of the expected value of V, and 
the standard deviation of N, are calculated. The expected value of N., N,, is 


given by 





(26) NV.) = = ap Pall), 


and the standard deviation of N;, oy,, is given by 
(27) on. = NE —(N,)’, 
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a 


From equation (8) it is seen that P,,(t) vanishes when t < nT, or when n > ¢/T. 
Thus the sums in (26) and (28) are, in fact, finite series although they are 
written as infinite series. If the last term of the series is of the form 0/0 (that 
is, £ = nT for some n) the term is to be taken equal to zero. This assumption 
makes NV, and N2 continuous functions of ¢ for t > 27. Now 


where 


(28) Ne = > 


n=0 








n 2, gttiy™ 
2 - 
(29) t—nT 


provided that t > nT. Since®P, vanishes for t < nT, the expression (29) may 
be substituted in (26) for all values of » without difficulty. The result is 


co co ao To 


co” Ze "at Py(l), 





~~ m+1 ’ 
m=0 t 


m=0 n=0 
or 
(30) v= Dae, 
where +! is defined by equation (7). Similarly, it is found that 
(31) i= ee 


An asymptotic expression for N, is obtained by substituting the result of 
equation (18) in equation (30). The result is 
(32) N.~ v—vT?/2t. 


Similarly, the asymptotic expression for V2 is 
(33) No ~ v4" (1t0T—r'T"). 


Finally, from (27), the asymptotic expansion for cy,’ is 


(34) one (1 + v7). 


IV. DISCUSSION 


According to equation (32), the expected, or mean, value of N, is vy, provided 
that the observation time, ¢, is large enough. This, of course, is to be expected 
since the purpose of calculating NV, is to obtain a close approximation to ». 

The value of the standard deviation of N, is more interesting. It might be 
thought that the standard deviation of the calculated value, N., is approxi- 
mated simply by (N,/t)!, because N, is the calculated rate of occurrence of a 
series of random events. However, equation (34) shows that this value must 
be multiplied by the correction factor (1+N,.T)! to take account of the addi- 
tional dispersion due to the dead time of the counter. Thus, a good approxi- 
mation to the standard deviation of N, is given by 
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" N.\ 

(35) oN, = (22) (1+N.T)' 
for large values of the observation time, t. 

It may also be remarked that, for most normal conditions, the correction 
term v?7?/2¢ in the expression (32) for NV, is much smaller than the standard 
deviation of N,. 

Finally, the method used here would be useful in other cases in which one 
wishes to determine the expected value of some function of the observed rate 
of registration, No. All that is necessary to apply the method is that the function 
in question should have a power series expansion whose radius of convergence 
is at least 7—!. 

APPENDIX 
It must still be shown that 
P = vt ee ; p” Ql =m vt te 
3) tin (2) | tn ft Ben (2) 
oe) 1+0T. eh O-\T7 
When this is shown, equation (16) follows. 
Equation (12) may be written in the form 





ce ver? on-m ye or m—1 —; 
(37) 3 Ln = Fy Anew) (m = 1,2,...). 


From equation (37) it follows that ™(t) satisfies the integral equation 
re: t m—1 — 
(38) n™(t) = i ge 2. mCz n (t—x)dx 
r k=0 


‘ — 
+ J ve" @—T n™(t—x)dx (t>T; m=1,2,...), 
= 


while 

(39) n™(t)'= 0 (0<t<T,m=1,2,...) 
and 

(40) n(t) = 1 (t > 0). 


First, it will be shown that n™ possesses a lower bound, given by (17), by 
the method used by Feller (1948) for the case m = 1. Consider the more general 
integral equation 


(41) A(t) = H(t)+ Taon A(t—x)dx (t> T). 
If 
(42) H(t) < {a > mC *(t—x)dx 


fort > T, and 
(43) A(t) <0 


for 0 << T, then consideration of the integral equation satisfied by 
n™(t) —A(t) shows that 
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(44) A(t) < n™(t) 
fort >0. Let 
_ Or ‘3 
(45) A(t) ” L 1+oT vr 
=0 (t < T). 
Then if A(t) is to satisfy (41), 
a pe=3 | (T’<#¢ 27), 
(46) H(t) " 
bon 7 —v(t—2T) av(t—2T) s a 
. [= -GpT J, ete ciate 


Since y*(t) = 0 fort < Tand k = 1,2,...,m—1, 
ot m—1 
(47) fer? & aCe mtt—x)dx = er 
YT 


for T <t<2T. Now, fort > T, 
eM) > 1+r(t-T), 


or 


* F4et—7)" 


and hence 





_ neem — _vt—T) | oe re 
7 * 1+v(t—T) * 1+v(t—T) (m = 1,2,...). 


For 7 <t < 2f, 
1+v(t—T) < 1+0T. 


Hence 


_ eticd ae _ 
(48) l—e > | ear (we = 1,2,...). 


Therefore inequality (42) holds for T < t < 27. 
Now assume that 


(49) n* (t) >) en oe (t> T; k= 0,1,...,m—1), 
and let 
2t—T m—1 ~ 
(50) G(t) = J oe?) = te |e | ee ie dx +- Si e77E-T gy 
T kno 1+vT 
Then 


t m—1 
(51) Git) < f ve" SF Cyn" (t—x)dx (t > 2T). 
T k=0 
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A change of variables in the integral shows that (50) may be written 
av(t—2T) 
ne ZB ht-8T) gr 
(52) GW) =e thes + J aha y; dy 


v(t—2T) vim 
ey 
= ——_-_—— d i 
J, (l-+0T)" 


Now, if it can be shown that G(t) — H(t) is positive or zero for ¢ > 27, then, 
because of inequality (51), inequality (42) holds. From (52) and (46), 


aise es vv(t—2T) m 
(53) G(t)—H(t) =e | 1c + i o(1+;25) iy| 


tT ; 
-[ 14oT (¢ > 27). 





s . ‘ a ee ee vl ast 
(54) G(2T)—H(2T) = l—-e 3 | >0 (ma = 1,2,...), 


because of inequality (48). Moreover, 


(55) 5 fe" (G@-HON 
“oe = m m—1 
= pot te- DP -bE-T )!"—m[ve-T)}"}. 


Now, if the first term in the braces on the right-hand side of equation (55) is 
expanded using the binomial theorem it is seen that the right-hand side is 
greater than or equal to zero for m > 1 and t > 2T. Thus, e”*?”[G(t) — H(#)] 
has a non-negative first derivative for ¢ > 27 and is itself non-negative for 
t = 27. Therefore G(t)— H(t) is non-negative for t > 27. This result, com- 
bined with inequality (51), shows that inequality (42) holds for t > 27. It was 
shown earlier that (42) holds for T < t < 27. Thus, inequality (44) holds, and 
therefore inequality (49) holds for k = m whenever it holds fork = 0,1,..., 
m—1. Since (49) holds for k = 0, it follows by induction that it holds for all R. 

It has now been established that A(t) is a lower bound to n™(t). (It seems 
possible that an upper bound could be determined in the same way but this 
question will not be dealt with here.) Since A(t) is a lower bound to ™(t), it 
can be shown that 


a Wee bes mTy™ m—1 
When n™—[vt/(1+vT)]™ satisfies this condition, a more general Tauberian 
theorem, given by Doetsch (1943), states that 


6 oy Ca) )- 8 SLR (Ga) J 


To show that the right-hand side of equation (57) may be replaced by 


ine Lro-(ia) |, 








it is necessary to modify a lemma given by Doetsch (1943, p. 209). 
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This lemma shows that when the integrand in (57) is monotonic increasing 
the indicated replacement may be made. The modification will take account 
of the fact that n™(t), rather than the complete integrand, is monotonic 
increasing. 

Since n™(t) is monotonic increasing, 





— “e ama 
(58) n™ (6t)(t—Ot) < J n”(s)ds < n™(t)(t—6t), 
6t 
for 0 < @< 1. Hence, it can be shown that 
here m ae m+1 \ 
(5¢ ieee m at. vt ) a. 
4 . )-\T4 07) Gtia—os 
£4 Leo Ge Io Loa I 
= wl ae 7) 1% 
a Ni _ut J. i of \ 
<r (4 (m+1)(1—6)S° 
Let 
; — m at _ sh vs. m 
(60) B= lim im J E s) =) Jas. 


Then, if f— © in the second part of the inequality (59), 


B(1—0") mv) = ( a) ee 

(61) m(1—6) samt \” (t) — lor (m+1)(1—6)$’ 
for0 <6@< 1. If 6— 1 in (61), the result is 

fc mint oof -(_} 
(62) B¢< lin inf € \” ™(t) idor) §° 
In a similar fashion it can be shown that 

> —(m— vI—m ay 
(63) E> lim sup f° \” ™(t)— = {- 

Inequalities (62) and (63) taken together show that 

. ~~) = =F is 

(64) lim ¢ \” (t)— = r f == B, 


Finally, equations (60), (64), and (57) show that equation (36) holds and hence 
that equation (16) holds. 
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PHOTOCONDUCTIVITY IN INDIUM-DOPED SiLICON' 


By J. S. BLAKEMORE 


ABSTRACT 


Impurity photoconductivity has been studied in silicon doped with indium, 
at various temperatures in the range 63-90°K. The long wave limit of spectral 
response and analysis of the bulk electrical properties both confirm an ionization 
energy of 0.16 ev. for the indium centers. An apparent variation of quantum 
efficiency in the range 5-8 yw is noted and discussed. The magnitude of the photo- 
conductive response suggests an effective recombination coefficient of about 
10-® cm./sec. at 90°K., rising steeply on further cooling. This demands a much 
more active recombination mechanism than present theories provide. 


1. INTRODUCTION 

The subject of photoconductivity by excitation of impurity levels in semi- 
conductors has received increased attention in recent years, starting from the 
work of Burstein et al. (1953) and of Rollin and Simmons (1953). The process 
is interesting both for the information it gives about the processes of carrier 
excitation and recombination, and for the possibility it presents of a sensitive 
infrared detector. 

The experiments mentioned above were carried out with silicon containing 
impurities of unknown type. In more recent investigations (as reviewed by 
Schultz and Morton 1955) material of known composition has in general been 
selected, since the impurity ionization energy depends on the chemical nature 
of the centers present. This ionization energy affects the long wave limit of 
photoconductivity and the magnitude of the cross-sections for excitation and 
recombination. 

A dependence of ionization energy on the nature of the center is not pre- 
dicted from the simple hydrogenic model for bound impurity states (e.g. 
Kittel 1953), according to which an electron or hole of effective mass m* in a 
semiconductor with dielectric constant K is bound to a singly charged impurity 
with an energy E = 13.6(m*/m)K~ electron volts. More detailed theories 
(Kittel and Mitchell 1954; Kohn and Luttinger 1955) consider cases where the 
tensor of effective mass does not reduce to a scalar, but the dependence of 
ionization energy on the nature of the impurity ion remains unexplained. 

For the elements of Group III of the periodic table as acceptor type im- 
purities in silicon, the ionization energies (Burton 1954) range from 0.045 ev. 
with boron up to 0.16 ev. with indium. With boron, aluminum, or gallium, a 
photoconductive response would be expected to persist for wavelengths up to 
about 20 yu, and the optimum temperatures for observing such effects would lie 
in the liquid hydrogen or liquid helium regions. With indium as the dominant 
impurity, photoconductivity is to be expected for wavelengths shorter than 
about 8 yu, and good conditions for observation of this should be developed in 
the liquid air temperature range. The latter subject was chosen for investi- 
gation. 
1Manuscript received April 23, 1956. 

Contribution from the Physics Department, University of British Columbia, Vancouver, 
British Columbia. 
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2. EXPERIMENTAL ARRANGEMENT 
(i) Necessary Conditions for Observation of Photoconductivity 


Since the acceptor levels provided by indium centers in silicon lie about 
0.16 ev. above the valence band, it is important that the photoconductive 
process should not be adversely affected by the presence of lower lying acceptor 
levels, arising from impurities such as boron. For if these lower levels remain 
partially ionized at the operating temperature, the rate at which free holes are 
thermally excited from them will greatly exceed rates of thermal or optical 
ionization for the more remote indium levels, and photoconductivity from the 
latter will be extremely small. The criterion is essentially that photoconduc- 
tivity may be expected if at low temperatures the Fermi level converges on the 
indium level, but not if the Fermi level tends towards any lower acceptor level. 

It is not necessary to rely on the improbable event of a crystal in which 
boron and other elements of similar behavior are completely absent, since 
such centers will be ineffective provided they are kept permanently ionized 
by the simultaneous presence of an equal density of donor centers, which may 
be deliberately introduced for this purpose. In practice exact cancellation is 
impossible and one must expect to use crystals in which there is a small surplus 
of donor centers over low lying acceptors. Electrons from the surplus donors 
ionize an equal number of indium levels. This raises the Fermi level, lowers 
the conductivity of the semiconductor, and affects photoconductivity to an 
extent dependent on the nature of the effective recombination mechanism. 

Crystals of indium-doped silicon were received from several sources for this 
investigation. Some of these revealed in their properties the presence of un- 
compensated boron and showed no photoconduction at 90° K. The most 
successful specimens came from crystals provided by Dr. R. Newman, in 
which some antimony had been deliberately included. These crystals were 
similar to those used by Newman (1955) in studies of this material. 

For observation of photoconductivity, the semiconductor should be main- 
tained at a temperature where the thermal rate of carrier excitation is not too 
large compared with optical generation for an obtainable illumination. With 
the materials and techniques employed here, the optical generation rate could 
at best be made some 107! of the thermal rate at 90° K., and smaller than this 
by a further factor of a thousand at the long wave limit. This serves to illus- 
trate the extent to which random fluctuations of the voltage developed across 
the specimen (i.e. noise) must be suppressed if meaningful results are to be 
obtained. 

(11) Experimental Techniques 

In our arrangement, the specimen holder was mounted inside a vertical 
glass tube, the lower end of which was sealed, whilst the upper end was fitted 
with a KRS-5 window. This specimen tube was surrounded with liquid oxygen 
or nitrogen in a glass Dewar, and the liquid refrigerant could be further cooled 
by forced evaporation. Thermal contact between liquid and specimen was 
maintained by filling the specimen tube with helium gas at a moderate pressure. 
The level of refrigerant was kept well above that of the specimen in its tube, 
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so that only a very small solid angle for the specimen was subtended by objects 
at room temperature. Thus, once equilibrium was reached, the specimen was 
very nearly at the temperature of the cooling agent. 

The radiation from a carbon arc was used in conjunction with a Perkin 
Elmer model 83 monochromator with NaCl prism to provide a monochromatic 
beam which was focused on the specimen, through the KRS-5 window. The 
radiation was mechanically interrupted (usually at 870 c.p.s.) and the a-c. 
photoconductive output at the chopping frequency measured with a tuned 
detector. 

The influence of current noise in the contacts was minimized by use of a 
four contact arrangement (Burgess 1955 a) and by employing essentially 
constant current conditions, current being derived from a d-c. source with a 
resistance much higher than that of the specimen itself. The output was 
obtained from the two potential probes, which themselves carried no current. 
These were connected to a preamplifier of very high input impedance, which 
led in turn to a General Radio wave analyzer. 

The design of specimens was also aimed at reduction of contact noise. The 
four contacts were applied by soldering with indium, using a method due to 
L. Young, details of which will be published elsewhere. In this way contacts 
were developed which were substantially ohmic and comparatively free from 
excess noise at the operating temperatures. The contact current density was 
minimized for a given bulk current density by etching the central section of a 
specimen down to a small cross-section, whilst the contacts were made on the 
end regions of much larger cross-section. 


3. DISCUSSION OF RESULTS 


(i) Spectral Response of Photoconductivity 

The general behavior of the spectral response of photoconductivity for 
indium-doped silicon is illustrated in Fig. 1 for a specimen operated at tempera- 
tures of 90° K. and 77.5° K. The ordinate in this figure is the responsivity R, 
defined as the fractional change in conductivity for unit energy density of 
illumination, R = (A¢/oo/), expressed in the units cm.?/watt. 

It will be seen that the intrinsic response for photon energies greater than 
1 ev. is much larger than the impurity response. Thus if any short wavelength 
light is scattered in the monochromator and emerges in the exit beam when 
measurements of the long wave response are being made, the accuracy of such 
results could be seriously affected. It was established that this did not occur 
for Our measurements, since the responsivity in the extrinsic range was not 
affected by interposing short wave absorbing filters of silicon or germanium 
at the monochromator entrance slit. 

The pronounced rise in responsivity for the extrinsic process on going to- 
wards smaller values of photon energy is primarily a consequence of the 
situation that each absorbed photon of energy greater than the threshold value 
will usually create only one free hole. Thus in Fig. 2, the extrinsic photo- 
conductivity per incident photon is given for the same data as that of Fig. 1. 
The equal-photon response still shows a dependence on photon energy, which 
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Fic. 1. Equal-energy responsivity R = (Ao/ool) cm.?/watt in the intrinsic and extrinsic 
response regions for indium-doped silicon (specimen 68). 


Equal-photon Response (Arbitrary Units) 
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Fic. 2. Extrinsic equal-photon response for specimen 68. (Thickness = 0.043 cm., (Na — Ni) 
= 5X10'§ cm.~) 
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derives from a variation of quantum efficiency and of absorption constant 
with the quality of the radiation. The quantum efficiency, Q, is the number of 
holes produced per absorbed photon ; which may be almost unity for most of 
the spectral range under consideration, but will fall for photon energies around 
the threshold value (=0.16 ev.), when other absorption processes become more 
important. We expect that the manner in which the quantum efficiency falls 
off at the extrinsic long wave limit will be broadly the same for all specimens 
of indium-doped silicon. Thus the general features of Fig. 2 agree with the 
equal-photon response of Newman (1955) taken at 77° K. However, a detailed 
examination must take account of the variation with hv of the proportion of 


incident photons which are absorbed. 


(11) Photoconduction per Photon Absorbed ; Quantum Efficiency 
Of an incident photon beam, a fraction r is directly reflected from the first 
f 
division of these depends on the thickness ¢ of the specimen, and on the absorp- 
tion coefficient a. For an incident beam of J watts/cm.?, the average volume 
rate of photon absorption throughout the semiconductor is 
(1) F = I. f(at)/1.6X10-'%(hv) photons/cm.’ sec. 


) 


surface of a specimen, and the remainder are transmitted or absorbed. The 


with hy expressed in electron volts. The function f(at), allowing for multiple 
reflections inside the specimen, is 
(2) OL) a ee cm?" 

Pict t [l—r-exp(—at)] 
This reduces to a when at < 1, whilst for strong absorption such that at > 1, 
it tends towards (1 —r)/t. 

These expressions enable us to relate the change in free hole concentration 
Ap for a given illumination with the quantum efficiency Q and the extrinsic 
hole lifetime 7. Clearly, 

(3) Ap = FQr cm. 


where 7 is the average lifetime for which a carrier is effective in conduction 
before suffering recombination. We can alternatively express Ap in terms of 
the density of holes in the unilluminated state, po, and the fractional increase 
in conductivity (Ag/oo) which results from an illumination; and hence in 
terms of the responsivity R. For 


(4) Ap = po(Aa/ao) = polR cm.=3 


where R is the equal energy responsivity as defined earlier. From (1), (3), and 
(4), the product of lifetime and quantum efficiency is 
(5 Or = Lox 10-9 poR(hv) flat) sec. 


In order to evaluate this product, @ and py must be known for the material 
with which the responsivity is measured. Since Newman (1955) has measured 
a as a function of fy for one impurity concentration, we may compute @ and 


f(at) for any sample, assuming @ to be proportional to the density of neutral 


acceptors. 
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By analysis of conductivity and Hall effect data, we may follow the variation 
of po with temperature from which impurity concentrations can be deduced. 
For if the semiconductor contains NV, indium levels per cm.* at an energy €a 
above the valence band, of which JN; are ionized at all temperatures with 
electrons from donor impurities, then fp» is given by the quadratic expression 


2 = j y Nie a \ y IV a 
(6) Po + poy it 8 exp( -¢4) "oe (V,—-N;) B exp( -.) = 0 


where 8 is the acceptor occupancy spin factor and N, = 2(2amkT7/h*) 


> 


«m*/m)*” is the effective density of states in the valence band. In p-type 
silicon, the mass parameter (m*/m) appears to be about 0.5 (Morin and Maita 
1954). 

The variation of free hole concentration py with reciprocal temperature is 
shown in Fig. 3 for the material which was used for our specimen 68, whose 


10'”, 


10'° 


2 4 6 8 10 12 
1oo0osT (%-') 


Fic. 3. Variation of free hole concentration with temperature for indium-doped silicon. 


response was illustrated in Figs. | and 2. This curve can be fitted on the 
assumptions that V, = 5.1X10'® indium atoms per cm.’, with V; = 8.510" 
cm.~* of these compensated by donors, and an acceptor ionization energy 
€a = 0.16(0) ev. 

Using the appropriate values of (Vz—N;,) and of po, the product Qr of 
Equation (5) can be calculated over the extrinsic spectrum for specimen 68, 
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using the observed responsivity and a reflection coefficient of 0.30.* The result 
is shown in Fig. 4, where Q7 is expressed as a function of wavelength rather 
than of photon energy so that the low energy region is displayed in the maxi- 
mum detail. This figure reveals a number of interesting features, regarding 
both spectral form and magnitude. 






QT (seconds) 
> 


Wavelength (microns) 


Fic. 4. Spectral variation of the product of quantum efficiency Q and effective hole lifetime 
7, for specimen 68. 


It seems probable that the lifetime of a hole will not depend to any marked 
extent on the energy of the photon which liberated it, since the excess kinetic 
energy is dissipated during the first few collisions with the lattice, in a time of 
the order of 10-'* sec. Thus the spectral dependence of Q7 may reasonably be 
ascribed to be that of Q itself. This appears to be constant in the wavelength 
range 1.8-5 u (hv = 0.25-0.7 ev.) and it is probable that the quantum efficiency 
is unity here, corresponding with lifetimes of 1.210~-® sec. at 90° K. and 
2010+" sec-atd1.0° i 

If this conclusion is drawn, then the slight rise in Q7 which occurs at 90° K. 
for wavelengths shorter than 1.8 yu implies a quantum efficiency here of more 
than unity. This could occur through impact ionization of another acceptor 
by a newly liberated hole, moving with velocity v = [2(hv—e,)/m*]! = 
7X107 cm./sec. for radiation of 1.3 4 wavelength, where the increase in Q is 
about 20%. This is consistent with a cross-section for impact ionization 
Simp = 6X10-'4 cm.2 (a not unreasonable value), if we set (Q—1) = 
vtoo imp(Na—N;), where to is the mean free time. 

The spontaneous impact ionization which occurs in germanium at very low 
temperatures above a critical field strength (Burstein et al. 1954; Ryder et al. 
1954) was not present in our measurements. This was expected, since with the 


*This follows from the expression r = (n—1)?/(n+1)? for reflection at normal incidence, 
where the refractive index n = 3.4. The influence of the extinction coefficient on reflection is 
completely negligible in the extrinsic spectral region. 
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field strengths employed, hole drift velocities were only about 10° cm./sec., 
much smaller than the mean thermal velocity in the non-degenerate hole gas, 
d = (2kT/m*)? = 7X10® cm./sec. at 90° K., which in turn is an order of 
magnitude smaller than the hole velocities which appear capable of promoting 
impact ionization. 

The slow fall in Q7 beyond 5 u suggests that the quantum efficiency begins 
to drop below unity, and is down to about 0.45 at 7.7 u (hv = 0.16 ev.), which 
marks the beginning of a much more rapid fall in efficiency. This behavior is 
seen most clearly at 90° K. and appears to be followed at 77.5° K., though 
the onset of noise prevented observations beyond 7.4 uw at this temperature. 
More limited measurements on other specimens at 90° K. and 63-64° K. also 
suggest the same type of behavior. 

It is not immediately obvious what other absorption mechanism is causing 
the fall in efficiency in the range 5-7.7 uw. Free carrier absorption can certainly 
be neglected, and there is no prominent lattice absorption band in this region. 
It might be argued that high energy photons may excite holes into the /} 
valence band (which is separated from the valence edge by the spin-orbit 
splitting), whereas lower energy photons may only create free holes in the 
ps2 bands which form the upper limit of the valence band, Tos, at k = (000) 
(Herman 1954). However, the small spin-orbit splitting of 0.035 ev. in silicon 
(Herman 1955) would appear to place the boundary of any such effect at 
about 6.4 u, whilst in fact a change in Qr is evident from 5.0 » onwards. 


(iii) Recombination of Carriers 

It seems likely that unit quantum efficiency is obtained in the spectral 
region 2—5 yw, and measurements of Q7 for this region are accordingly taken 
as indicating the effective lifetime for holes. (This is the lifetime for extrinsic 
capture of holes. The time scale for intrinsic recombination of holes and elec- 
trons will be quite different.) The extrinsic recombination rate 1/7 is propor- 
tional to the density of ionized acceptors (= N; for weak illumination at low 
temperatures where po < N;), the constant of proportionality being the 
recombination coefficient B, cm.*/sec. This quantity is determined by the 
cross-section of each ionized acceptor for capture of a free hole, ¢, cm.*, and 
by the mean thermal velocity of free holes, i cm./sec. These relationships are: 
(7) 1/r = B,N; = o,0N; sec.—. 


B, and o, are characteristic of the recombination process and not functions 
of impurity densities. Thus these are the quantities which may be compared 
with theory. Values of B, deduced from the magnitude of Q7 for three speci- 
mens are shown in Fig. 5 for the various operating temperatures. (Note that 
a recombination coefficient of 10~* cm.*/sec. corresponds with a capture cross- 
section of some 10~" cm.?) In view of the large temperature dependence, a 
semilogarithmic representation in terms of 1/7 seemed appropriate, and the 
line for B, « exp(0.055/kT) appears reasonably consistent with the general 
trend. The internal accuracy of the various points in Fig. 5 is not very high, as 
is not surprising in view of the number of separately measured quantities which 
are involved in converting measurements of signal voltage for each specimen 
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Fic. 5. Recombination coefficient B, for three specimens, expressed as a function of 
temperature. 


into the equivalent values of B,. Nevertheless, the results are sufficiently clear 
to demonstrate their incompatibility with the extant theories of recombination 
(Sclar and Burstein 1955 ; Gummel and Lax 1955) both in respect of magnitude 
and of temperature dependence. 

These theories are all based on the hydrogenic model for bound impurity 
states, which is patently inapplicable for silicon containing indium, where the 
ionization energy of 0.16 ev. is some three times larger than that expected 
from the model. Thus no better than order of magnitude agreement between 
theory and experiment for B, was to be expected in any case. In fact, the 
measure of disagreement is more than three orders of magnitude. 

The two processes of optical recombination and impact recombination dis- 
cussed by Sclar and Burstein yield recombination coefficients of 410-% 
cm.*/sec. and 107'! cm./sec. respectively for our case at 90° K. The tempera- 
ture dependence is small for the first model whilst in the second case B, is 
proportional to the free hole concentration, and will accordingly be even smaller 
at lower temperatures. It is clear that neither of these mechanisms can play a 
dominant part in the present circumstances. 

A rather larger recombination rate of some 3X10-'° cm.*/sec. would be 
predicted for the thermal recombination model of Gummel and Lax if one- 
phonon processes were permissible. But no phonons in silicon have energies 
greater than about 0.106 ev. (Hsieh 1954), and for transitions of free holes to 
the ground state of indium impurities, we must have recourse to multi-phonon 
processes, with correspondingly smaller probability, particularly at low 
temperatures. 





oa 
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This wide lack of agreement between experimental measures of recombina- 
tion in germanium and silicon and the above theoretical models was noted by 
Sclar and Burstein in connection with results for impurities of low activation 
energy at 4.2° K. They felt that the principal cause of the discrepancy lay, not 
in the inadequacy of the hydrogenic model, but in the presence of another 
unknown mechanism for recombination which has not yet been considered. 
The present results strongly endorse this conclusion, and point to one char- 
acteristic of this process in particular; that it seems to show a very large 
negative temperature coefficient. 

It is interesting to note that an energy of 0.055 ev. arises in the interpre- 
tation of Fig. 5. This is about the energy we should expect for the compensated 
shallow acceptor levels above the valence band. Thus it is possible that tempo- 
rary deionization of some of these states may be significant in the recombina- 
tion process. Alternatively, it is possible that free holes first relax to an excited 
bound orbit around an indium ion, subsequently progressing to the ground 
state. These carriers would be lost to the photoconductive process as soon as 
they relapsed into bound.states of high quantum number, and observations at 
low intensity would yield no information about multi-stage recombination. 
To investigate such possibilities, it would be necessary to study non-linearity 
of photoconductivity at very high light intensities, when intermediate states 
might become saturated. As an alternative, carriers in various orbits might 
be observable using paramagnetic resonance. 


(iv) Detectivity 

The ease with which low intensity infrared radiation may be studied is 
determined by the extent to which noise in the specimen and detection system 
is excluded. The use of the techniques outlined in Section 2 served to minimize 
this interference, but even so the noise was considerably higher than the level 
set by thermal and shot noise (Burgess 19555). For the best specimen studied 
(No. 68), the excess noise factor was some 15 db. at 870 c.p.s. when operated 
at 90° K., and more at lower temperatures. 

The very rapid increase in noise voltage when specimens were cooled below 
90° K. caused the detectivity to be poorer at the lower temperatures, even 
though the equal power responsivity was enhanced. It is possible that 
optimum conditions for study would have been found at a temperature slightly 
higher than 90° K., so that both signal and specimen noise would be reduced, 
the latter by the largest factor. The permissible rise in temperature could only 
be a small one, since noise in the detection system would soon take over as the 
limiting factor. 

The detectivity was, of course, a function of wavelength, reaching a maxi- 
mum in the neighborhood of 5 uv. For operation at 90° K., the detectivity of 
specimen 68 (area 0.035 cm.’, thickness 0.043 cm.) at this wavelength was 
1.0109 watts! for a one cycle bandwidth at 870 c.p.s. This is roughly 
comparable with the detectivity quoted by Rollin and Simmons (1953) for 
silicon at 14° K. containing acceptors of smaller activation energy (possibly 


boron). 
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NEUTRON YIELDS FROM ACTINIDE-BERYLLIUM ALLOYS! 


By O. J. C. RUNNALLS AND R. R. BOUCHER 


ABSTRACT 


The neutron yields from the reaction Be%(a, z)C™ have been measured for 
several actinide—beryllium alloys. The values for the beryllium thick target 
neutron yield, derived from the measurements on alloys containing Pu®*, Am*!, 
and Cm*?, are 60+2, 74+2, and 112+3 neutrons per million alphas, respectively, 
when based on the Chalk River standard Ra”®-a@-Be source. The dependence 
of the thick target yield, mmax, on the alpha-particle energy, E, is described by 
the empirical formula mmax = 0.152 E%-® neutrons per million alphas. The yields 
for Ra-a-Be and Ac”?-a-Be, calculated from the latter expression, are in reason- 

‘ able agreement with the best observed yields from powder-compacted sources. 
For Po-a-Be, the calculated thick target yield of 67 neutrons per million alphas 
lies within the limits of error of the yield observed by several earlier workers. 


INTRODUCTION 


Neutron sources have been prepared in past years by mixing a salt of an 
alpha-emitting radioactive isotope with a light element. The highest neutron 
yield is obtained when beryllium is employed as the target nucleus from the 
reaction Be®(a, ~)C! (Anderson 1948). The usual source of alpha-particles 
has been either Ra??® or Po?!®. Recently, however, several alpha-emitting 
elements of the actinide series, notably plutonium, americium, and curium, 
have become available as by-products from the operation of nuclear reactors. 

The actinide elements are relatively non-volatile when compared with 
radium or polonium and all probably form alloys with beryllium. Thus, the 
fabrication of alloy neutron sources rather than pressed mechanical mixtures 
would appear feasible, and of some advantage in reproducibly approaching 
the maximum neutron yield. 

A knowledge of the variation of the neutron yield with the beryllium: 
actinide atom ratio in such alloys should yield, on extrapolation, a value of the 
thick target neutron yield for the Be®(a, 2)C!? reaction. The energy dependence 
of the thick target yield should also be measurable, using alloys containing 
single alpha-emitters of differing energies; e.g. Pu?*®, Am**!, and Cm?#?. 


EXPERIMENTAL 


(a) Preparation of Alloys 
The available thermodynamic data (Brewer et al. 1949, 1950; Sense e¢ al. 
1954) indicate that the actinide halides may be reduced by beryllium to metal 
if the reactants are heated in vacuum; e.g., for the reaction 
2PuF 44) +3Be(.) > 2Puc:) +3BeF x) (1) 


the calculated pressure of gaseous BeF, over the system at 1400°K. is approxi- 

mately 20 mm. Hg. If excess Be is added initially, so that the|product is a 

beryllium—plutonium alloy, the equilibrium will be shifted further to the right 

in equation (1) because of the reduction in the activity of liquid plutonium. 
1Manuscript received April 4, 1956. 
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Alloys of plutonium, americium, curium, and actinium with beryllium were 
prepared in BeO crucibles by the reduction of the actinide trifluoride with 
powdered Be metal, 99.8% pure, in vacuum at 1100°-1200°C. At these tem- 
peratures the beryllium fluoride product readily distilled, leaving a fluoride- 
free alloy. The Be-rich alloys were melted by heating to 1375°C. 

The plutonium trifluoride, 99.5% pure, was prepared by fluorinating 
plutonium dioxide with a mixture of dry HF and Hz at 650°C. in a Pt tray. 
The specific activity of the plutonium used was 1.58 X10° d.p.m./mgm., i.e. 
the Pu?4° content was 4.8%. 

The Pu content of each alloy, 100-300 mgm., was determined from the weight 
of the PuF; added, on the assumption that the reduction yield was 99%. 
The latter figure was obtained after several alloys had been dissolved and 
assayed for plutonium by alpha-counting. The Be content was determined by 
difference, after the alloy was weighed. 

Two americium-beryllium alloys were prepared by the Be reduction of 
30 mgm. amounts of AmF3. The alpha-radioactive purity of the americium 
was >99.99%. The total inactive impurities were <0.5%. The Am**! content 
of one alloy was estimated from the starting weight of AmF;, with a suitable 
correction for the Am lost by distillation during the alloying process, as 
described earlier (Runnalls and Boucher 1955). The Am content of the second 
alloy was estimated by comparing its gamma-intensity with that of the first 
alloy. 

Two curium-beryllium alloys were prepared by the Be reduction of 20-30 
ugm. amounts of CmFs, carried on 30 mgm. of PuF3. The radioactive purity 
of the Cm**? was 99.5%, before the plutonium carrier was added. The Cm 
contents of the alloys were determined by a-counting, after dissolving the 
small ingots. A suitable correction, about 2%, was made for the Pu added. 

An actinium-beryllium alloy was prepared from 13.3 mc. of Ac*?? which 
was in radioactive equilibrium with its daughters. The Ac was mixed with 
50 mgm. of carrier uranium, and the mixture fluorinated at 650°C. in a Pt 
tray. The resulting fluoride was heated with powdered Be metal in vacuum 
for 30 min. at 1100°C., before raising the temperature to 1350°C. for 15 min. 
The product was a single shiny ingot which had shrunk away from the BeO 
crucible on solidification. Because of the appreciable gamma-activity from the 
above quantity of Ac it was possible to monitor for losses which occurred 
during evaporation, fluorination, alloying, etc. An estimated 98%, 13.0 mc., 
of the original Ac and 49 mgm. U were present in the final 455 mgm. alloy 
ingot. The Be:Ac+U atom ratio was, therefore, 245. 


(b) Neutron Counting Methods 

The alloy sources were sealed in aluminum cans 3 cm. long and 2 cm. 
in diameter before they were removed from the glove-box system. The neutron 
yields from the canned alloys were determined by a direct comparison of their 
neutron output with that from standard Ac-a-Be and Pu-a-Be sources using 
the simple paraffin block counter described by Hanna and Runnalls (1956). 
Although such a counter will be sensitive to neutron energy, it was assumed 
that the neutron energy spectra for the different actinide-alpha-beryllium 
sources would be similar. The Ac and Pu sources had been calibrated against 
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the Chalk River Ra-a-Be standard and a National Bureau of Standards 
Ra-y-Be source by activation of a MnSQ, solution, as described in the follow- 
ing paper (Hanna and Runnalls 1956). 

The chemical reaction between the single alpha-emitting fluorides and bery!l- 
lium could be followed by positioning the neutron counter near the reduction 
furnace, since the neutron yield increased as the alloy formed. The change 
in neutron yield on the solidification of a liquid alloy could also be observed, 
and related to the absolute yield by a later comparison of the cooled, solid 
alloy with a standard source. 

RESULTS 
(a) Plutonium-Beryllium Alloys 

The neutron yields from 11 Pu—Be solid alloy sources are plotted as circles 
in Fig. 1. Several of the alloys had shown a higher neutron output when 
molten, as indicated by the squares in Fig. 1. The dashed curve, representing 
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Fic. 1. The variation of the neutron yield with atom ratio for plutonium—beryllium alloys. 
the theoretical yield, has been constructed from the following expression, 
suggested by Baerg (unpublished), relating atom ratio and neutron yield, 
i.e. 


n/Nmax = N . Spe/(N . SpetSru)s 


where 
n/Nmax = fraction of maximum theoretical neutron yield, 
N = Be:Pu atom ratio, 
Spe = Be atomic stopping power = 0.63, 
Spy = Pu atomic stopping power = 4.6. 
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The values quoted for the atomic stopping powers relative to air were ob- 
tained by extrapolating known data for alpha-particles of about 6 Mev. 
energy (Livingston and Bethe 1937). Reasonable agreement between the 
theoretical neutron yields and the highest observed values is obtained when 
Nmax, the maximum theoretical neutron yield, is set equal to 60.0 neutrons per 
million alphas, as Fig. 1 illustrates. 

The neutron yields were appreciably affected by the crystallite size and the 
extent of agglomeration of the only intermetallic phase, PuBeis, which formed 
in the Pu-Be alloys. A plot of the neutron emission rate during the beryllium 
reduction of two PuF; charges is shown in Fig. 2. At 1125°C. the alloy forma- 
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Fic. 2. Measurements of the neutron emission rate during the conversion of two plutonium 
trifluoride — beryllium mixtures to sintered, and subsequently molten, plutonium-beryllium 
alloys. The effect of solidification on the neutron emission rates from molten alloys A and B 
is also shown. 


tion rate was rapid, as indicated by a rapid increase in the neutron output. 
The chemical reduction appeared complete after 25 min., since the neutron 
output became constant. The product at this stage consisted of a sintered 
mixture of PuBei; aggregates in a Be matrix. On raising the temperature to 
1375°C., where the two alloys A and B were molten, the neutron yield in- 
creased abruptly. When alloy A solidified no change in the neutron output was. 
observed. When alloy B solidified, however, the neutron yield decreased 
markedly as illustrated in Fig. 2. The Be:Pu atom ratios in alloys A and B 
were 265:1 and 285:1, respectively. Metallographic examination showed that 
the PuBe,; crystal size was ~2 yu in alloy A and ~ 15 u in alloy B (see Fig. 3). 
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Fic. 3. Polished sections of plutonium—beryllium alloys A, 90.9 wt.% Be, and B, 91.5 
wt.% Be, as-cooled, after etching for 30 sec. with 1% HF, showing dark PuBe1; crystals in 
the light Be matrix. 

Fic. 4. Polished section of an 89.2 wt.% beryllium-uranium alloy containing 13.0 mc. of 
Ac®7, as-cooled, after etching for 30 sec. with 1% HF, and showing the dark UBeis phase in 
the light Be matrix. 


The range of an alpha-particle of initial energy 5.14 Mev. in a PuBe,; 
crystal may be calculated assuming that the range is inversely proportional 
to the stopping power of the target, i.e. 


Reuse; = Rate . Nair/(N puseis . S pues)» 


where 
Rair = range of 5.14 Mev. alpha in air = 3.65 cm. (Livingston and 
Bethe 1937), 


Nair = atoms per cm.? of air = 5.110'® at 15°C. and 76 cm. Hg, 
N puse,; = Molecules of PuBeis; per cm.? = 7.4X10?!, calculated from the 
X-ray density of 4.35 gm. cm.—* (Runnalls 1956), 
Spuse,;; = molecular stopping power of PuBeis; = 4.6+13 X0.63 = 12.8. 
Thus the calculated range of the plutonium alpha in PuBe,; is 
Reupeis ~ 20 u. 
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It follows that the larger the PuBe; crystals in Be-rich alloys are, the closer 
the neutron yield will approach that of a 13:1 atom mixture of Be:Pu, as an 
increasing fraction of the alphas are trapped in PuBe,s. 

The measured neutron yield for PuBeis; was 38.7 neutrons per million 
alphas or 63% of that obtained for alloys with Be:Pu ratios of ~250. If one 
assumes that the specific activity of the Pu used by Stewart (1955) was 
1.4108 d.p.m./mgm., the neutron yield from her PuBe3 source was ~40 
neutrons per million alphas. 

(6) Americium-—Beryllium Alloys 

An Am*4!-Be alloy with a Be-Am atom ratio of 263:1 did not decrease in 
neutron yield during solidification from the liquid, as already described by 
Runnalls and Boucher (1955). Presumably, the AmBe,; crystals were com- 
parable in size to the PuBe,; crystals in Pu-Be source A referred to previously, 
i.e. ~2 yu. The measured neutron yield was 71.7 neutrons per million alphas. 
The neutron yield from a second source with a Be:Am atom ratio of 14:1 was 
48.5 neutrons per million alphas. The calculated maximum theoretical yields 
were 73.7 and 74 neutrons per million alphas respectively, for a homogeneous 
mixture of atoms. 


(c) Curium-Beryllium Alloys 


The neutron yields from the two curium-beryllium alloys are indicated in 
Table I. Since neither molten alloy showed any apparent decrease in neutron 














TABLE I 
NEUTRON YIELDS FROM TWO CURIUM-BERYLLIUM ALLOYS 
Be: Cm+Pu Total Cm** Total Cm neutron Neutrons 
Source atom a-activity a-activity neutron yield per million 
ratio (d.p.s.) (d.p.s.) yields (n. /sec.) Cm alphas 
(n./sec.) 
1 217:1 2.00X10® 1.95109 2.17XK108 2.14X105 110 
2 240:1 5.32X10® 5.20K10® 5.61105 5.55105 107 


yield during solidification, it is assumed that the carrier PuBe,; crystals were 
~2 wu diameter, or less. The calculated maximum theoretical yields, assuming 
a curium atomic stopping power of 4.7, are 114 and 110 neutrons per million 
alphas for sources 1 and 2, respectively. 


(d) Actinium—Beryllium Alloy 

During the beryllium reduction of the AcF;-UF, mixture, an appreciable 
fraction of the Ra?** and short-lived daughters distilled from the system, 
resulting in a decrease in neutron yield. Thus, the increase in neutron emission 
rate accompanying the alloy formation was masked. When the melt solidified, 
the neutron output increased slightly as the short-lived daughters from the 
Ra*?s remaining in the alloy grew in and were trapped in the solid. As a result, 
the effect of UBei; crystallization on the neutron yield was not directly 
observable. 
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A sample of the distillate, recovered from the surface of the glass bulb 
containing the furnace assembly, was examined on an alpha-pulse analyzer. 
The only activities observed were those due to Ra?** and its daughters, i.e. 
no Th??? distilled. This was confirmed by the later observation that the 
neutron yield from the alloy increased with the half-life of Ra??*, 11.2 days. 
After 90 days, the neutron yield had increased to a constant value of 2.85 x 105 
n./sec., a value determined from the activation of a MnSO, solution (Hanna 
and Runnalls 1956). 

A pressed Ac?*?-a-Be source, prepared by the Commercial Products Division 
of Atomic Energy of Canada Ltd. from a part of the same batch of Ac solution 
described above, was compared in both neutron and gamma output with the 
alloy. The pressed source had been prepared by evaporating a slurried mixture 
of Ac(OH); and —325 mesh Be powder to dryness, followed by compacting 
the powder under pressure in a steel die. The source was contained in a Monel 
capsule with a silver-soldered lid. The Be:Ac atom ratio was approximately 
65,000:1. 

The ratio in gamma output between the pressed source and the alloy source 
at equilibrium was 0.995; i.e. the Ac??? content of the pressed source was 
12.9 mc. compared to 13.0 me. in the alloy source. The neutron yield from the 
pressed source, however, was significantly higher, 3.25 10° n./sec. 

The lower yield, 2.85 10° n./sec., in the alloy source was attributed to the 
presence in the alloy of UBe,;-carrier* crystals appreciably larger than 2 uz. 
The latter supposition was confirmed when the alloy was metallographically 
examined (see Fig. 4). The UBe:; crystals containing the Ac were, in fact, 
50-100u in diameter. 


DISCUSSION 


The thick target neutron yields, 60, 74, and 112 neutrons per million alphas, 
for Pu?%*, Am*4!, and Cm?*?, respectively, are plotted in Fig. 5, as a function 
of their respective energies, 5.14, 5.47, and 6.10 Mev. The presence of 4.8% 
Pu?!° in the Pu does not affect the result since the average energy of the Pu?4° 
alpha is 5.15 Mev. (Asaro and Perlman 1952). Earlier measurements on thick 
target neutron yields by Segré and Wiegand (1944) and Roberts (1944), 
using Po?!® sources, are included in Fig. 5, together with the yield calculated 
by Anderson (1948) from Halpern’s data (1948). 

The possible errors on the mean thick target yields for Pu, Am, and Cm are 
estimated to be +3% owing to uncertainties in the determination of the 
actinide content of the alloys and in the comparative neutron measurements. 
They are subject to a further deviation of +5% because of the uncertainty 
in the absolute neutron yield of the Chalk River Ra-a-Be standard (Hanna 
and Runnalls 1956). The values quoted by Segré and Wiegand (1944) were 
believed absolute within 30%. The greatest possible error in Roberts’ value 
(1944) was stated to be due to a +20% uncertainty in the neutron emission 

*Since actinides and lanthanides such as Th, U, Np, Pu, Am, La, and Ce all form an XBei;- 
type compound with beryllium, presumably an isostructural lattice exists in Ac—Be alloys. 


Hence in the UBe:; crystals shown in Fig. 4 the Ac atoms probably occupy normal U positions 
in the UBe,; lattice. 
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Fic. 5. Thick target neutron yields on beryllium for Pu®*, Am™!, and Cm”? plotted against 
their respective alpha-particle energies. Thick target measurements by earlier workers using 
Po"? are included. 
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rate from the standard sources used. The possible errors in Halpern's work 
(1948) were not stated. 

An empirical expression relating the thick target neutron yield on beryllium, 
Nmax, to the alpha-particle energy, E, may be derived from the data for Pu**® 
and Cm?%?, as follows: 

Mmax = 0.152 E**> neutrons per million alphas. 


Using the above function, shown as a solid line in Fig. 5, the calculated max 
for Am**! is 75 neutrons per million alphas. The calculated thick target yields 
on beryllium from the Ra?** and Ac??? chains are listed in Tables II and III. 
They are in reasonable agreement with the best observed yields from powder- 
compacted sources. 


TABLE II 


Ra-a-Be THICK TARGET NEUTRON YIELDS 
CALCULATED FROM THE EMPIRICAL FORMULA 
Nmax = 0.152 E*® NEUTRONS PER MILLION ALPHAS 








a-Emitter Energy, E, neutron yield, 
Mev. Nmaxy 
n./108 a's 
Ra%6 4.777 46 
Em? 5.486 76 
Po?!8 5.998 105 
Po?! 7.680 258 
Po?!0 5.298 67 


For Ra-a-Be, the calculated yield varies from 485 to 552 neutrons per 
million alphas, or 1.79X107 to 2.04107 n./sec./gm. Ra, depending on the 
Po?! content of the radium. The maximum observed yield obtained from 
pressed Ra-a-Be sources has been quoted as 1.7X10' n./sec./gm. Ra by 
Anderson (1948), for freshly-separated radium, when no Po?!® is present. A 
somewhat higher yield of 1.9107 n./sec./gm. Ra is, however, regularly 
attained by the Commercial Products Division of Atomic Energy of Canada 
Ltd. in the preparation of small, pressed Ra-a-Be sources using radium 
containing not more than 50% of the equilibrium amount of Po?!® (Anstee 
unpublished). 

The thick target yield of 67 neutrons per million alphas, calculated for 
Po-a-Be, lies within the limits of error of the measurements of Roberts (1944) 
and Segré and Wiegand (1944), where the observed yields were 80+16 and 
73422 neutrons per million alphas respectively. The value of 50 neutrons 
per million alphas obtained by Anderson (1948) from Halpern’s data (1948) 
would appear somewhat low, when compared with the latter observations. 

The calculated thick target yield for Ac?*’-a-Be at equilibrium is 709 neu- 
trons per million alphas, or 3.38 10° n./sec. for 12.9 mc. of Ac???, only 4% 
higher than the observed yield of 3.25X10° n./sec. in the pressed source 
described previously. The latter high yield may be explained by assuming 
that the actinium hydroxide has deposited as a thin film over the available 
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TABLE III 


Ac”?-@-Be THICK TARGET NEUTRON YIELDS 
CALCULATED FROM THE EMPIRICAL FORMULA 
Mmax = 0.152 E*® NEUTRONS PER MILLION ALPHAS 


Thick target 





a-Emitter Energy, E, neutron yield, 
Mev. Nines, 
n./108 a's 
Th? 5.90 99 
Ra223 5.65 84 
Em?!9 6.71 158 
0215 7.365 222 
Bit! 6.56 146 


Be surface. If the average Be particle size for the latter source was 20 u and a 
density of 5 was assumed for the 12.9 mc. of actinium hydroxide, the calculated 
thickness of a uniform film on all the exposed surface is only 6 A. 

When molten Pu-Be alloys with a Be:Pu atom ratio near 265:1 solidify it is 
not fully understood why the resulting PuBe.; crystals are an order of magni- 
tude smaller than those forming in alloys with a higher Be:Pu ratio. It is 
suspected, however, that alloys in the former composition range undergo a 
monotectic transformation during cooling, whereas alloys richer in Be precipi- 
tate primary crystals of PuBei; during the whole period of solidification. 
A careful determination of the Be-rich end of the Pu-Be alloy equilibrium 
diagram would lead to a better understanding of the latter phenomenon. 
Unfortunately, in the only reported investigation of the Pu-Be system (Kono- 
beevsky 1955) little work on Be-rich alloys was carried out. 
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MEASUREMENT OF THE NEUTRON YIELDS FROM 
ACTINIDE-BERYLLIUM ALLOYS! 


By G. C. HANNA AND O. J. C. RUNNALLS 


ABSTRACT 


The neutron yields of a series of actinide—beryllium sources have been meas- 
ured. Some of the sources were compared with primary neutron standards in a 
manganese sulphate bath. The experimental procedure is described in detail and 
corrections for the loss of neutrons from the bath are discussed. These measure- 
ments include a comparison of the Chalk River Ra-a-Be and the National 
Bureau of Standards Ra-y-Be standard sources. The majority of the actinide- 
beryllium alloy sources were then compared with the standardized sources in 
simple boron trifluoride counter systems. 


INTRODUCTION 


Two different aspects of the neutron yields from actinide—beryllium alloys 
are of interest, the change in yield with a particle energy and composition of 
alloy, and its absolute value. The absolute scale of the measurements has been 
established by standardizing some of the sources against the Chalk River 
Ra-a-Be and the National Bureau’ of Standards Ra-y-Be standards in a 
manganese sulphate bath. A simple counting system, consisting of a boron 
trifluoride counter embedded in a paraffin wax block, has been used to compare 
neutron yields. Some comparisons have also been made using a long counter. 

1. THE MANGANESE SULPHATE BATH 
l(a) General 

There are various methods available for the accurate comparison of neutron 
sources having different energy spectra and certain of these are more or less 
directly capable of measuring the absolute value of the neutron yield (see 
for example Feld 1953). It was not however the intention of the present work 
to attempt an independent absolute determination, but rather to rely, at 
least for the present, on the standard sources available to us. 

The most accurate comparisons are probably those made using the well- 
known “‘integration’’ technique, in which the resonance, or thermal, neutron 
density distribution in a paraffin oil bath surrounding the source is explored 
using thin foils or miniature counters. The density integral for an infinite 
system is then computed. Water is often used instead of oil, which results in a 
greater energy dependence, fast neutron reactions ((, a) and (m, p)) being 
more important in oxygen than carbon with neutron spectra of the Ra-a-Be 
type. 

The measurement of the neutron density, using a finite number of detectors, 
and the subsequent numerical integration, is, however, a rather tedious 
process. In the manganese-bath method (Anderson, Fermi, and Szilard 1939) 
a continuous detector is used, and the integration is performed mechanically 


1Manuscript received April 18, 1956. 
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by stirring the manganese solution after activation; the specific activity of the 
solution is then proportional to the source strength. The method suffers, 
however, from the inherent drawback of an aqueous moderator noted 
previously. 

In its most sensitive form the manganese-bath method employs the Szilard— 
Chalmers reaction to concentrate the activity induced in a permanganate 
solution in a filterable precipitate of manganese dioxide, which is then 
counted as a compact source (Dodson, Goldblatt, and Sullivan 1946; Erber, 
Rieder, and Broda 1950; Edge and Shire, 1955). This technique is rather 
laborious, however, and in the present work the activity of the manganese 
solution was measured by immersing a scintillation counter in it. Descriptions 
of manganese-bath measurements of this type have been given by Taschek 
and Turner (1945), Almqvist (1950), McRonald (1951), and Littler et al. 
(1951). 

A disadvantage of the manganese-bath method compared with the density 
distribution method is that a larger volume of moderator is required. The 
neutron density distribution need be measured in only one direction (and hence 
the source may be placed non-centrally), and it can be extrapolated to correct 
for the loss of neutrons. In planning a manganese-bath experiment it is 
necessary to compromise between loss of neutrons, which cannot be estimated 
very accurately, and sensitivity, which is approximately inversely proportional 
to the total volume of the bath. 


1(b) Description of the Apparatus 

A cylindrical stainless steel tank 24 in. in diameter and 30 in. deep, with a 
9 in. diameter hole in the top, was filled to a depth of 24 in. with a concentrated 
solution of manganese sulphate. This solution, of specific gravity 1.37, con- 
tained 430 gm. of the anhydrous salt per liter. It is calculated that 51% of the 
neutrons were captured in the manganese. 

The tank was set up in a nominally inactive location well separated from 
other ‘‘counting’’ areas, but, unfortunately, not sufficiently well separated 
or shielded from various fluctuating sources of radiation. A considerable 
improvement in the background level was obtained by building a lead wall 
2 in. thick and 40 in. high around the tank. 

Each neutron source, delivered to the building by truck when required, 
was removed from its protective lead castle behind a lead wall, where it was 
placed in a suitable aluminum wire cradle attached to a long string, and 
whence it was swung rapidly over and into the center of the manganese bath. 

The scintillation counter consisted of a 2 in. diameter by 1 in. thick sodium- 
iodide crystal, a 5819 photomultiplier tube, and a cathode follower housed in a 
water-tight aluminum tube 3 in. in diameter and 2 ft. long. The output 
pulses from this counter were fed through long cables to a fast linear amplifier 
and discriminator, and counted on a Marconi scaler equipped with an auto- 
matic timing mechanism. 

The counter was used in three locations: (i) in the manganese bath, where 
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it was clamped in a standard position with the crystal at the center of the 
bath, (ii) in the background monitoring enclosure (described below), and 
(iii) to check the counter sensitivity, on a lead-shielded test bench at a standard 
distance from a Co® source. 

Preliminary experiments on the effect of bias on the sensitivity to background 
on the one hand, and the radiations from the activated manganese solution 
on the other, indicated that the best discriminator bias would be at about 
100 kev. It appeared later that a smaller bias would have been a wiser choice, 
for, although it would have increased the ratio of the background counting 
rate to that of the activated manganese without a useful increase in absolute 
sensitivity, it would have reduced the dependence of sensitivity on the over-all 
gain of the system, which was much less stable than expected. 

Apart from the actual counting of manganese activities, to be described 
below, two subsidiary investigations had to be made. First, the dependence 
of counter sensitivity on the gain was determined for both the Co*® source 
and the manganese solution. By deliberately altering the photomultiplier H.T. 
voltage it was found that the change of gain producing a 1% increase in 
Co® counting rate produced a 2.2% increase in that of the manganese solution. 

Second, it was necessary to be able to estimate the background during a 
measurement. It has already been noted that the background varied errati- 
cally in time; it was also found to vary considerably over quite modest dis- 
tances. To minimize the possible effect of the time variations differing between 
different points, the background monitoring station, a small, lead-shielded 
enclosure, was placed close to the manganese bath, outside the lead shielding. 
Ideally the counter, when placed in this enclosure, would have the same shield- 
ing as it had in the manganese bath. No attempt was made to fulfill this 
accurately; however, it was established experimentally that, within the accu- 
racy of the measurements, a unique relationship existed between the background 
measured in the shielded enclosure and that in the manganese bath. This 
relationship was checked at intervals during the main series of measurements 
at times when the manganese activity had decayed to a very low (and 
calculable) level. 


l(c) Experimental Procedure 

The activation measurements were made in the following way. In the case 
of the weaker sources, the neutron source was left suspended in the bath 
overnight to bring the manganese activity within about 2% of saturation; 
the stronger sources were exposed for appropriately shorter times. After 
removal of the source, a powerful electric stirrer was lowered into the bath, 
and stirring continued for 20 min., by which time the source was safely out 
of the building. 

Counting was then started beginning with a measurement of the Co® 
counting rate (three one-minute counts, ~60,000 c.p.m. in a background of 
~1000 c.p.m.). This was followed immediately by a measurement of the 
manganese activity (usually 10 one-minute counts), bracketed by two meas- 
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urements of the rate in the background monitoring station. This in turn was 
followed by another Co® check, and the routine continued until (usually) 
four manganese determinations had been made. The complete series of 
measurements lasted about two hours. 

The background in the bath was about 650 c.p.m.; the manganese activity 
was some 2800 c.p.m. for the smallest source (6.6 X10‘ n /sec.), and between 
10! and 2.510‘ c.p.m. for the others. 

Each measurement of the manganese activity was corrected for the back- 
ground in the bath, which was calculated from the mean of the two measure- 
ments in the background monitoring station taken immediately before and 
after the manganese measurement. The net counting rate due to manganese 
was then normalized to a standard Co® counting rate; the mean value of the 
two Co® measurements made before and after the manganese measurement 
was used. Since it could not be assumed that the counter sensitivity was 
changing steadily, the difference between the two Co® measurements repre- 
sented an uncertainty, often the major one, in the normalized rate. Finally 
this rate was corrected to infinite irradiation time and zero subsequent 
decay. 

1(d) Results 


Usually four measurements were available for each source. In all cases they 
were within +1% of the mean value, except for the smallest source, where the 
background correction was much more serious; the spread in this case was 
+2%. 

The series of measurements on five sources lasted one week, during which 
time the sensitivity of the counter to the manganese activity had fallen by 8%. 
As a check on the normalization procedure the first source was then remeasured. 
The value obtained differed by only 0.6% from the original one. 

In the case of the Ra-y-Be source virtually all the neutrons are slowed 
down to low velocities and die in the tank. The a-Be sources however give 
much higher energy neutrons and corrections are required for escape from the 
tank and for loss of neutrons in high energy reactions. These corrections, 
which are difficult to estimate accurately, are discussed in the Appendix. 
Approximately 4% of the neutrons escape from the tank and 3% are lost in 
fast neutron ((”, a) and (n, p)) reactions in oxygen and sulphur. 

The absorption of thermal neutrons in the material of the sources was 
estimated to be negligible in all cases except for the Ra-y-Be source where the 
correction was 0.4% and for the Pu*%a-Be source where a correction of 
—0.7% for the increased emission due to fission was applied. 

The final results are given in Table I. All strengths are given in terms of 
that of the Chalk River standard 200 mgm. Ra-a-Be source which was 
2.92 10° n./sec. at the time of the measurement. This value was obtained 
by correcting for the RaD that had grown in since the standardization of the 
source at the Argonne National Laboratory (June 21, 1950), where it was 
calibrated, with an accuracy of +2%, against A.N.L. source No. 38, which is. 
known absolutely to +7% (Seid! and Harris 1947). 
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TABLE I 
Source Neutrons/sec. Per cent error* 
Chalk River 200 mgm. Ra-a-Be standard (2.92108) —_ 
Pu Be No. 4 (625 mgm. Pu) 6.47 X10* +2 
Ac Be A (13.0 mc. Ac”) alloy source 2.85105 +1 
Ac Be B (12.9 mc. Ac”?) pressed source 3.25105 +1 
N.B.S. 1 gm. Ra-y-Be standard No. 2 1.26 10° +3 


*Not including the uncertainty of the Ra-a-Be standard source. 
The N.B.S. standard source has been described by Curtiss and Carson (1949). Further 
details of the plutonium and actinium sources are given by Runnalls and Boucher (1956). 


1(e) Discussion of Results 

The strength of the N.B.S. source No. 1 is 1.255108 n./sec. (De Juren and 
Chin 1955); N.B.S. No. 2 was originally (1/0.9871) times stronger (Curtiss 
and Carson 1949), but, having sprung a leak and been repaired since that time, 
it is now only 0.948 times as strong, i.e. 1.19010° n./sec. (Caswell 1956). 
The standard error quoted is +2%. A comparison between N.B.S. No. 2 and 
A.N.L. No. 38 has been made in a subcritical pile (De Juren, Padgett, and 
Curtiss 1955) which indicates that the A.N.L. scale of source strengths is 
3+1% higher than that of the N.B.S. Hughes (1954), considering all the 
available international intercomparison data, concludes that A.N.L./N.B.S. = 
1.05; the compilation of De Troyer and Tavernier (1954) suggests a value 
close to 1.03. It is to be expected then that the present measurement of the 
N.B.S. source, in terms of a source calibrated to +2% at A.N.L., would give 
a value that is about 4+3% higher than the N.B.S. calibration. The result, 
in fact, is 643%. 

It appears then that the Chalk River standard is about 5% too high with 
respect to the N.B.S. standard, which lies very close to the international ‘‘best 
value’ (Hughes 1954). In spite of this all source strength measurements 
have been expressed in terms of the Chalk River standard. This scale appears 
to be in reasonable agreement with that in use at the Commercial Products 
Division of the Atomic Energy of Canada where the strength of the Ac Be 
source B was independently measured as 3.22105 n./sec. (Anstee 1955). 


2. LONG COUNTER MEASUREMENTS 


A long counter was used to compare five small plutonium beryllium neutron 
sources with the larger sources that had already been standardized in the 
manganese bath. The long counter, developed by Hanson and McKibben 
(1947), consists of a cylindrical BF;-filled proportional counter embedded 
coaxially in a cylindrical assembly of paraffin wax with dimensions and general 
arrangement chosen empirically to achieve a substantially constant sensi- 
tivity over a wide range of neutron energy. 

The sources were mounted on a light framework some two feet from the 
front of the counter about five feet above the floor well away from the walls 
of the room and with all other equipment well behind the counter. The efh- 
ciency of the counter, 2 ft. from the source, was approximately 5X 10-°. 

The angular distribution of neutrons was not investigated in complete 
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detail but the observed asymmetries were small in all cases and would not 
have introduced errors of more than +2%. The results of the measurements 
of the small sources are not of interest here, but the results with the standard- 
ized sources are given in Table II, which shows the ratio of source strength 
measured with the long counter to that measured with the manganese bath 
(on both occasions relative to the Chalk River standard). The results are gen- 
erally satisfactory, but it will be noted that the long counter is more sensitive 
to the lower energy Ra-y-Be neutrons than to the a-Be sources by ~9+4%. 














TABLE II 
< Strength in long counter Strength in block counter 
Source 3 
Strength in Mn bath Strength in Mn bath 
Chalk River Standard (1.000) — 
Pu Be No. 4 1.01+.03 (1.000) 
Ac Be A 1.02+ .02; 1.01+.03 
Ac Be B 1.03+ .02; — 
N.B.S. Ra-y-Be standard 1.11+.04 —_ 


3. THE PARAFFIN BLOCK COUNTER 

Most of the routine source strength comparisons were made using a simpler 
and more sensitive counting system, consisting of a BF;-filled proportional 
counter embedded in a paraffin wax block 8 in. by 8 in. by 18 in. The counter 
was 2.5 in. in diameter and had an active length of 12 in. It was filled to a 
pressure of 60 cm. Hg with enriched BF; (96% B!°). The sources were placed 
at a distance of 20 in. from the long face of the block, giving an efficiency of 
about 107%. 

This system will be sensitive to neutron energy, but it was felt that the 
effects of energy differences between Pu®*, Am?4!, and Cm? would be small. 
Measurements were always made using both Pu Be No. 4 and Ac Be A stand- 
ards. Table IJ shows that the block counter gave a ratio for these two sources 
in good agreement with that obtained in the manganese bath. 

The angular distribution of neutrons was expected to be rather similar in 
all the sources and it was therefore ignored, all measurements being made 
with the axes of the sources parallel to the counter axis. 

ACKNOWLEDGMENTS 

The authors are indebted to Dr. W. G. Cross for the loan of his long counter 
system, to Mr. C. G. Tario for much help with the manganese bath measure- 
ments and to a large number of their colleagues for many valuable 


discussions. 


APPENDIX: THE EFFICIENCY OF THE MANGANESE BATH 
In comparing neutron sources of different energy spectra in a manganese 
bath, it is necessary to correct for the differential loss of neutrons (cf. Sections 
l(a), 1(d)). There are minor differences between the different (a, m) type 
sources, but these have been ignored, and an attempt has been made to 
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correct only for the much more pronounced difference between Ra-a-Be and 
Ra-y-Be sources. In the latter case, the neutron energies are low enough 
for the loss to be negligible, and an estimate is therefore required of the loss of 
Ra-a-Be neutrons in high energy neutron reactions and by escape from the 
bath. 

For Ra-a-Be neutrons, which have a maximum energy of about 12 Mev. 
and a mean energy of about 5 Mev. (Teucher 1949), the significant fast neutron 
reactions are O'§(n, a), S®(n, a), and S®(n, p). It would not be easy to calculate 
the frequency of these reactions even if the cross sections were known over the 
required energy range; in fact, data exist only below about 4 Mev. (Hughes 
and Harvey 1955). De Troyer and Tavernier (1954) have measured directly 
the loss of Ra-a-Be neutrons due to fast neutron reactions in oxygen as 2.6+ 
0.3%. For S®, a plausible extrapolation to higher energies of the cross section 
data plus a very crude calculation suggests that, in the bath used, where the 
S/H atom ratio was 0.03, about 1% of the neutrons are lost in (m, a) and (n, p) 
reactions. 

The escape of neutrons from the bath could, unlike the loss due to fast 
neutron reactions, be made vanishingly small by using a very large bath, 
but at the cost of a reduced sensitivity. The bath used in the present measure- 
ments, which represents a compromise between escape and sensitivity, was a 
cylindrical volume of nearly saturated manganese sulphate solution 24 in. in 
diameter, 24 in. deep. 

The escape of neutrons cannot be calculated from the age theory of slowing 
down. This theory is (subject to other limitations that are also invalidated 
here, see Marshak (1947)) applicable only up to a distance from the source 
of about L,?/X,, which, for Ra-a-Be neutrons in water, is about 7 ¢m. (Zs, 
the slowing down length, and \,, the scattering mean free path, are both 
about 7 cm.). Beyond this distance, neutrons that have travelled directly 
from the source without making a collision become increasingly important, 
and the slow neutron distribution, Gaussian on the age theory, becomes 
exponential in character. Marshak (1947), Verde and Wick (1947), Wick 
(1949), and Holte (1954) have developed improved theories of slowing down, 
but the numerical evaluation of their expressions would be extremely laborious. 

It is possible, however, to derive a sufficiently accurate estimate of the 
escape from experimental data. There have been several determinations 
of the thermal neutron density p surrounding a Ra-a-Be source in water 
(Tittle ef a/. 1949). Beyond about 15 cm. from the source, pr? is a quite 
accurately exponential function of the distance 7; thus the fraction of neutrons 
dying beyond r is A exp(—r/Le). The measurements of Rush (1948) give 
A = 1.91 and Le = 9.47% cm. 

The escape of neutrons from a bare sphere of radius 7 will be greater than 
this by a boundary correction factor, B, say. The estimation of B is difficult. 
If the density of thermal neutron sources remained unaffected up to 7, the 
increase in escape would be due solely to thermal diffusion, giving for water 
B, = 1.2». This source density is, however, less than in the infinite bath 
because of the escape of those fast neutrons that would have been reflected 
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back from beyond 7; this effect is probably comparable with that of thermal 
diffusion, giving B = 1.4 (water) and B’ = 1.3 (manganese sulphate). It is 
interesting to compare this estimate for a 23 cm. radius sphere of water, 
viz. 1.4X1.91 Xexp(—2.427) = 0.236, with the experimentally measured 
value of 0.22+0.05 obtained in a 40 gm./liter solution of potassium per- 
manganate by Edge and Shire (1955). 

For a nearly saturated manganese sulphate solution, the escape is expected 
to be B’A’exp(—r/Le’), where B’ = 1.3 and A’ = A = 1.9. The solution 
used in the present work (density 1.37 gm./ml.) contains some 5% fewer 
hydrogen atoms per unit volume than ordinary water, which would make 
Le’ = 1.05 Le; on the other hand, inelastic scattering in manganese and 
sulphur will counteract this effect. A satisfactory calculation of the inelastic 
scattering appears to be impracticable. Very roughly, however, since Le is 
determined essentially by the scattering cross section of hydrogen for the 
primary neutrons (roughly 1.5 barns), then, with Mn/H = S/H = 0.03, a 
5% decrease in Lp would result from inelastic cross sections of mw(R+4A)?, 
1.5 barns for Mn, 1 barn for S. Thus the two effects compensate and Lp’ = Lr 
within two or three per cent, substantially independent of the manganese 
sulphate concentration. 

Rumsey and Ward (1948), using strong manganese sulphate solutions 
(300 gm./liter), have measured the relative efficiency of cylindrical baths of 
different sizes, the largest being 24 in. deep, 24 in. diameter. The expected 
escape from a cylindrical volume may be obtained by numerical integration 
of the exponential formula, with a significant uncertainty in the effective 
value of r only near the corner of the cylinder, where the escape is least 
(maximum r). Rumsey and Ward’s data for 24 in., 20 in., 16 in., and 12 in. 
cylinders are fitted well by the integration of the expression B’A’exp(—r Lp’) 
with B’ = 1.3, A’ = 1.9, and Le’ = 8.43 cm. = 0.89 La. It seems unlikely 
that inelastic scattering is three times as important as the crude estimate 
made above, as it would have to be to explain such a low value of Lz’. Possibly 
A’ is less than A because of a marked difference in the neutron distribution 
near the source, and B’ may have been overestimated. With Lp’ set equal to 
Lr, the expected value, good agreement with the experimental data is obtained 
with B’A’ = 1.9. 

There are then three estimates of the escape from a bare sphere of radius r 
filled with a 300-400 gm. /liter solution of manganese sulphate: 

(i) 1.3 A exp(—r/Lp) 

(11) A exp(—r/Lp) with A = 1:0. 2, = 9.47,-cm:, 

(iii) 1.3 A exp(—r/0.89 L»)| 


where (1) is calculated from the observed neutron distribution in water, and 
(ii) and (ii) are obtained from fitting the Rumsey-Ward data to exponential 
formulae. 

The percentage of neutrons escaping from a 24 in. diameter, 24 in. deep 
cylinder of manganese sulphate solution calculated from these expressions is: 
(i) 7.8%, (ii) 6.0%, (iii) 5.38%. 
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The present bath was surrounded by a lead wall 2 in. thick, 40 in. high, 
open top and bottom, which, it is calculated, would have reflected back into the 
bath 50% of the neutrons hitting it. The net escape is calculated to be: 
(i) 4.7%, (ii) 3.6%, (iii) 3.1%. 

Allowing for uncertainties in the numerical integration and the calculation 
of the effect of the lead wall, and compromising between the different estimates 
of escape, a value of 4+1% for the net escape is obtained. 

To this must be added the loss due to high energy neutron reactions, and 
therefore the total correction to be applied when a Ra-a-Be source is compared 


with a low energy one amounts to 7+2%. 
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LATITUDE EFFECT OF THE COSMIC RAY NUCLEON AND MESON 
COMPONENTS AT SEA LEVEL FROM THE ARCTIC TO THE 
ANTARCTIC' 


By D. C. Rose, K. B. Fenton, J. KatzMAN, AND J. A. SIMPSON 


ABSTRACT 


Results are presented of cosmic ray measurements taken at sea level during 
1954-55 from the Arctic to the Antarctic. The equipment consisted of a neutron 
monitor and a meson telescope. Latitude effects of 1.77 for the nucleonic com- 
t and 1.15 for the meson component were measured. The longitude effect 
- equator was much less than expected on the basis of the geomagnetic 
ric dipole and the longitude effect at intermediate northern latitudes 
hows that the 1 tude of the effective eccentric dipole is considerably west of 
that of the geomagnetic eccentric dipole. In a previous paper by the same 

ithors, the positions of the equatorial minima were combined with other 
published cosmic ray measurements to calculate a new cosmic ray geomagnetic 
equator. In this paper new coordinates are derived on the assumption that these 

| coordinates apply to a new eccentric dipole, and, therefore, that the 
| coordinates may be extended to high latitudes. When the complete 
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res plotted on these coordinates, it is found that an eccentric dipole 
re of the earth's magnetic field is inconsistent with the combined 
bservations at all latitudes. 


I. INTRODUCTION 


For over twenty years there has been evidence that the observed distribution 
of cosmic ray intensities is not in agreement with simple Stérmer theory for 
the motion of cosmic ray particles in the dipole field of the earth. 

The early meson component observations were difficult to interpret, how- 
ever, because of both the small latitude effect and the difficulty of correcting 
for the atmospheric temperature effects. In this paper the properties of the 
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ce level 


ititude and longitude effects have been studied by measuring the 


leonic component intensity which is free of atmospheric temperature 


effects and displays a very large latitude effect. The positions of the lower 
energy limits (the ‘knees’ of the latitude curves) for primary particle detec- 


tion have been measured in both hemispheres. Measurements were made at 


the same time with a meson counter telescope in order to compare these new 
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During 1954-55 there were two expeditions extending from 89°N. to 75°S. 


yeomagnetic latitude. Through the co-operation of the Royal Canadian Navy 
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d the United State Navy the equipment was carried on board the ice 
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breaker H.M.C.S. Labrador during its circumnavigation of the North American 
continent (23 July to 20 November 1954) and on board the icebreaker U.S.S. 
{tka « ts Antarctic expedition from Boston (1 December 1954 to 138 April 
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CRUISES OF LABRADOR AND ATKA TAKING COSMIC RAY MEASUREMENTS 


Fic. 1. Voyages of H.M.C.S. Labrador and U.S.S. Atka. The northern voyage of the 
Labrador (solid line) was from Halifax, north through the North West Passage, south along the 
Pacific Coast of North America, through the Panama Canal, then east to the Island of Grenada, 
and straight north to Halifax. The southern voyage of the Atka (broken line) was from Boston 
directly to the Panama Canal, then to New Zealand, south to Antarctica, east around the 
coast of Antarctica, northwest to Buenos Aires and Rio de Janeiro, and back to Boston. 


The dates covered by these observations include the end of the solar mini- 
mum and the beginning of increased solar activity which was responsible 
for modest cosmic radiation intensity-time variations during January and 
February 1955. Consequently intensity variations with time have not ser- 
iously interfered with the interpretation of the data. 

During the entire period of observations, fixed station monitors were in 
operation at Ottawa and at Climax, Colorado. 

Il. THE NUCLEON COMPONENT 


A. Measurements and Corrections 

The equipment carried on the two expeditions is described by Rose and 
Katzman (1956) in their presentation of the results of the northern cruise. 
(This paper will be referred to as Reterence A.) The equipment consisted 
of a neutron monitor of the type developed by Simpson but with some modi- 
fications in design for use on shipboard. The neutron counting rate at high 
latitudes was about 10,000 counts per hour and two separate recording chan- 
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nels were used so that a continuous check on its operation could be maintained. 
The equipment was similarly mounted in both ships, there being a } in. 
steel deck overhead and about 4 in. of thermal insulation. In the Labrador, 
the mast structure was directly overhead but this appeared to have a negligible 
effect as would be expected since it covered a relatively small solid angle. 
The same equipment was used for both expeditions. 

The data have been corrected for atmospheric pressure changes and world- 
wide cosmic ray intensity variations. 

For pressure correcting the results of both cruises, a pressure coefficient of 
—0.735% per mb. was used, this being the best value selected from many 
analyses of data from fixed stations carrried out at Ottawa and Chicago. 
Pressure coefficients were worked out using the present data for latitudes 
above the knees in both the Arctic and Antarctic, the values being —0.74+ 
0.01% per mb. in the north and —0.67+40.05% per mb. in the south. The 
results selected for the southern analysis were obtained during a period when 
there were minor intensity variations of the type usually associated with solar 
activity. A low pressure coefficient and of low accuracy is quite common 
under such circumstances. 
The corrections for world-wide intensity variations were made in the 


following way. The variations at Ottawa, observed with an identical neutron 
monitor, were reduced to percentages of the mean intensity for the whole 
period (after pressure corrections had been made). The data were plotted and 
slightly smoothed, and proportional corrections were then made on the 
pressure-corrected results taken at sea. Although it is now established that 
the amplitude of world-wide variations is dependent upon particle rigidity, 
the most important intensity variations occurred when the Atka was at high 
southern geomagnetic latitudes (January and February 1955). Consequently, 
using data from Ottawa (geomagnetic latitude 56.8°N.) for corrections should 


not introduce appreciable error. The corrections were small (being less than 
3°7) and have no effect on the over-all results of the latitude variation but are 


of some importance in considering both the normalization of the measurements 
aboard the two ships and differences in intensity between the Arctic 

and Antarctic well above the latitude knee. 
The results, after these corrections, are shown in Fig. 2 and a complete 
table of the data is included as Appendix I. Twelve-hour averages were taken 
nst the mean geomagnetic latitude during the 12 hours. Where 


; 1 j 
and niotred aga 
and P10 i agai 


the ships were stationary or the position was changing very slowly, particu- 
larly in the Arctic and Antarctic, longer averaging periods were used in order 


to avoid obscuring the details of some parts of the figures with too many points. 
The accuracies are discussed in Reference A. While the counting rate would 
indicate a standard deviation of 0.4% for the 12 hr. points, it is estimated 
that this should be raised to 0.4% because of other inaccuracies. 
No normalization has been carried out in Kig. 2. Though the same equip- 


ment was used and similarly mounted, small differences due to effects of 
surrounding parts of the ship might be expected. ‘two methods of normaliza- 


tion are possible. One is to use the data taken at Balboa where both ships 


(about a month apart) for two or three days. The other is to assume 
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that the mean intensities north of the northern knee and south of the southern 
knee should be the same. As plotted, the mean intensity in the north, using all 
the data above the knee, is 0.68% higher than that in the south. If the data 
had been normalized using the measurements at Balboa the northern intensity 
would be 3.59% higher than that in the south. These differences, even if real, 
do not in any way influence the conclusions in this paper. 

Fig. 3 shows the data around the knee in the southern hemisphere. Here 
corrections for the world-wide variations have been made using data taken 


160 






-counts per hour + 64 


Neutron intensity corrected for 
primary intensity variations 
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Fic. 3. The high latitude knees of the nucleon component latitude effect in the southern 
hemisphere. 
from the neutron monitor at Climax. A comparison of the intensity variations 
at Ottawa, Climax, and Chicago during January and February, 1955, showed 
them to be similar. In Fig. 3 some of the data taken while sailing eastward 
around Antarctica have been omitted. The figure is presented to show in 
more detail the irregular variations found after making the best corrections 
and the difficulty in locating very accurately the latitude of the knees. 


B. Discussion of Results 

Since the latitude effect and its variation with longitude are of great impor- 
tance in determining the effect of the earth’s magnetic field on cosmic rays, 
these will be discussed in some detail. The most striking feature is the mag- 
nitude of the effect for the nucleon component as measured by such a neutron 
monitor. Previous sea level measurements on the hard component show a 
latitude effect of about 1.1 (ratio of the intensity at high latitudes to that at 
the equator). The present data gave a latitude effect of 1.77. A large effect for 
the nucleon component was of course expected from the measurements of 
Simpson et al. (1951, 1953), taken at airplane altitudes. 

Since the curves in the north and south hemispheres are flat above approxi- 
mately 55° for the sea level measurements of the nucleon component, this 
method of detection is sensitive to primaries of rigidity as low as 1.6 Bv. 
This estimate is based on conventional geomagnetic coordinates and vertical 
critical rigidities. 
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The locations of the minima for the two latitude curves at 100°W. and 30°W. 
longitudes do not coincide with the geomagnetic equator. This problem has 
been studied in some detail elsewhere (Simpson ef al. 1956) where it is shown 
that these data and all other available cosmic ray data for intensity minima 
in equatorial regions may be fitted with a change in geomagnetic coordinates 
but without changing appreciably the tilt of the geomagnetic dipole. Later in 
this paper the effect of extrapolating these revised geomagnetic equator 
coordinates to high latitudes will be investigated. 

For any two latitude curves having different intensities at the equatorial 
minimum, a crossover point in both the northern and southern hemispheres 
would be expected where the intensity is essentially independent of longitude. 
This arises from the quadrupole contributions of the geomagnetic field, 
equivalent to an off-center dipole, and has been computed in detail by Jory 
(1956). In Fig. 4, two idealized latitude curves are drawn for different geo- 
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Fic. 4. Sketch of predicted form of latitude curves for different geomagnetic longitudes 
in an eccentric dipole field. 
magnetic longitudes showing the crossover near 35°N. and S. A comparison 
of this figure with Fig. 2 reveals that the behavior of the experimental data 
in the southern hemisphere is in qualitative agreement with the theory, but 
the latitude curves in the northern hemisphere display an inverted longitude 
effect without the crossover expected on the assumption that the earth’s mag- 
netic field is represented by an off-center dipole. 

In Fig. 5 the neutron intensity as a function of geographic longitude has 
been plotted for a number of different north geomagnetic latitudes including 
~35° where the crossover latitude region has a negligible longitude effect. It is 
clear from Fig. 5, however, that the longitude effect is not only large but passes 
through a maximum at ~70°W. longitude. Assuming that the correct coordin- 
ates were used, this change of sign for the slope of the longitude curve represents 
the longitude through which passes the meridian plane containing both the 
center of the earth and the magnetic center of the dipole (with the magnetic 
center farther from the surface). From geomagnetic theory, this plane has 
been located at ~20°W. longitude, but it has been suspected for many years 
(Hoerlin 1936) that the position of this plane must be appreciably west of 
20°W. longitude for the off-center dipole effective for cosmic ray particles. 
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Fic. 5. The longitude effect for several northern geomagnetic latitudes. The data used for 
these curves are obtained from Fig. 2 and Appendix I. 


The difference in cosmic ray intensity for the minima at 30°W. and 100°W. 
(Fig. 2) is ~2%, whereas if the longitude effect corresponded to the geo- 
magnetic coordinates the value should have been ~8%; this reduced longitude 
effect at the equator probably arises from the westward shift of the eccentric 
dipole as indicated by the data in Fig. 5. 

C. Assumed Geomagnetic Coordinates 

One possible way of removing some of the disagreements between the results 
in Fig. 2 and the expected result for an eccentric dipole field would be to try 
a new set of coordinates. Simpson et al. (1956) have shown that a cosmic ray 
geomagnetic equatorial plane appreciably different from that for the con- 
ventional geomagnetic dipole field gives a better representation of both these 
and other published data for the positions of the equatorial minima. Sup- 
porting evidence for this equatorial plane comes from measurements of the 
azimuthal asymmetry at high altitudes near the equator (Danielson et al. 1956). 
Furthermore, measurements of primary a-particle rigidities and fluxes at high 
altitudes by the nuclear emulsion technique give some confirmation of the view 
that new coordinates are needed up to geomagnetic latitudes of 30° or 40°N. 
We shall, therefore, consider a world-wide set of coordinates based on the 
position of the cosmic ray geomagnetic equatorial plane. 

The normal to the new equatorial plane intersects the surface of the earth 
in the north at latitude 80.2°N. and longitude 113.2°W. Using this axis, a 
new set of cosmic ray geomagnetic coordinates has been computed* to see if 
the present measurements lead to a better picture on the basis of an eccentric 
dipole field. The results are shown in Fig. 6, and the new coordinates are 
included in Appendix I. 

The minima are now close to the equator (the best fit, using all available 


data, did not quite agree with the Atlantic equator crossing) but the longitude 

* Note.—The magnetic coordinates were computed by the University of Toronto Computing 
Centre on the digital computer Ferut. A program for this machine has been prepared by Dr 
J. Hart in such a way that any axis could be used and either geographic coordinates inserted 


or a table derived for conversion of coordinates over the surface of the earth 
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effect at latitudes up to the knee is much larger than in Fig. 2 where conven- 
tional geomagnetic coordinates are used. Such a large longitude effect may not 
in itself be impossible but a comparison with Fig. 4 shows that the relative 
positions of the knees in the north and south are still not in agreement with a 
simple eccentric dipole field. If one assumes that the southern hemisphere 
data present the correct picture, the longitude effect in the north is in the 
wrong sense. Crossover points occur at about 20°S. and 10°N. which could at 
best be said to be only in qualitative agreement with Fig. 4. The situation 
between 10°N. and 30°N. is complicated by the fact that the route of the 
Labrador and the south-bound part of the U.S.S. Atka route were through 
the Panama area at these latitudes and the longitude was changing rapidly. 

The selection of the new cosmic ray geomagnetic coordinates has therefore 
still left the following difficulties: (a2) it has accentuated the longitude effect 
and in so doing has (0) accentuated the anomaly in the northern hemisphere. 


III. THE MESON COMPONENT 
The data obtained with the triple coincidence telescope of counting rate 
14,000 counts per hour are shown in Fig. 7 where conventional geomagnetic 
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Fic. 7. Variations of the meson component intensity with geomagnetic latitude. 


coordinates are used. The results have been pressure corrected using a pressure 
coefficient of —0.18% per mb. On the southern cruise, 24-hr. averages were 
used instead of 12-hr. averages. 

The maximum latitude effect is about 1.15 and since this is small and tem- 
perature corrections could not be made (upper atmosphere temperature data 
were not available during the northern voyage), the analysis of the meson data 
in relation to the geomagnetic field will not give nearly as satisfactory results 
as in the case of the nucleon component. 

There are a number of features worth mentioning, however. The most 
important has already been included in the paper by Simpson ef al. (1956) 
on the equator crossings, that is, the positions of the equatorial minima for 
the meson component and for the nucleon component agree. This agreement 


shows that the positions of the equatorial minima obtained by other observers 
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using meson recording equipment may be used with more confidence than was 
possible previously. However, the relative intensities at the Atlantic and 
Pacific equator crossings are in the opposite sense from those of the nucleon 
component. The surface temperature differences were in the right direction to 
indicate that this may be a temperature effect. 

The large difference between the North Atlantic and North Pacific intensities 
is discussed in detail in Reference A where it is shown that it is likely to be 
largely, though probably not totally, due to temperature. 

The latitude knees are not at all definite; in fact the meson data show no 
evidence of a sharp knee. The cutoff producing the knee would be expected, 
of course, to depend to a very large extent on the atmosphere. 


IV. CONCLUSIONS 


Even though the above discussion is based upon data derived from a small 
portion of the earth, it is already clear that a simple eccentric dipole magnetic 
field does not adequately describe the outer magnetic field distribution with 
which the incoming cosmic ray particles interact. Indeed, the evidence for a 
change of geomagnetic coordinates in the equatorial region already reported 
(Simpson et al. 1956) provides strong support for the view that the description 
of the outer magnetic field distribution is different from the description derived 
from surface magnetic field measurements. It is now known that the earth 
and its magnetic field are immersed in a highly conducting interplanetary 
medium, and it is probable that the equatorial effects, as well as the results 
reported here, arise from the interaction of the inclined rotating geomagnetic 
field with this ionized medium. It is possible that quadrupole and higher order 
terms of appreciable magnitude must be introduced to describe adequately 
the geomagnetic field which produces the observed distribution of cosmic ray 
interactions at the earth. 
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APPENDIX I 


Since the data taken during the two voyages to the Arctic and Antarctic may be of use to 
others to compare with past or future measurements they are presented herewith in detail. 
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COSMIC RAY DATA FROM Labrador AND Atka CRUISES 




















(A) Labrador 
Corrected 
Date Time, Pressure, _ results, Lat. 
U.T mb. counts/hour 
+64 

July 23 0-12 1002.5 : 44.5°) 
(1954) 12-24 1005.3 5.1 44.6 
24 0-12 1008 .3 5.8 46.3 
12-24 1013.3 5.3 48.7 

25 0-12 1019.5 3 51.2 
12-24 1022.0 3.7 53.5 

26 0-12 1021.8 3.7 56.3 
12-24 1021.2 5.4 59.8 

27 0-12 1021.3 164.0 61.3 
12-24 1024.6 166.7 64.2 

28 0-12 1026.6 167.0 66.9 
12-24 1026.7 167.0 69.5 

29 0-12 1026.0 166.2 q2.2 
12-24 1024.4 166.1 73.9 

30 0-12 1024.8 167.2 74.3 
12-24 1020.9 168.0 74.4 

31 0-12 1014.5 168.8 74.5 
12-24 1014.1 165.8 74.5 

Aug. 1 0-12 1015.4 168.0 74.5 
12-24 1016.5 166.6 74.5 

2 0-12 1017.4 166.5 74.4 
12-24 1015.8 166.0 74.4 

3 0-12 1014.3 166.6 74.3 
12-24 1016.6 167.4 73.8 

4 0-12 1017.0 169.1 74.2 
12-24 1017.5 166.8 74.5 

5 0-12 1019.1 167.3 74.7 
12-24 1020.8 167.0 75.6 

6 0-12 1022.9 166.0 75.9 
12-24 1021.3 165.0 75.6 

7 0-12 1022.8 167.8 76.4 
12-24 1024.1 164.7 7.1 

8 0-12 1022.9 165.6 78.7 
12-24 1023.0 166.1 78.9 

9 0-12 1025.2 167.1 78.9 
12-24 1024.9 167.8 78.9 

10 0-12 1020.9 167.7 78.9 
12-24 1020.5 168.3 78.9 

11 0-12 1024.7 166.8 78.9 
12-24 1024.9 167.8 78.9 

12 0-12 1023.7 166.8 78.8 
12-24 1023.8 165.8 78.2 

13 0-12 1019.9 167.3 76.9 
2-24 1019.3 165.4 76.1 

14 0-12 1015.6 167.0 75.4 
12-24 1012.6 167.8 74.5 

15 0-12 1012.4 168.3 74.3 
12-24 1009.5 168.6 74.6 

16 0-12 1010.5 168.5 74.6 
12-24 1010.0 168.6 74.6 

17 0-12 1008.6 168.0 74.6 
12-24 1008 . 2 167.5 74.5 

18 0-12 1012.2 166.7 74.0 
12-24 1013.6 167.8 74.6 

19 0-12 1012.9 167.7 74.5 
12-24 1011.9 166.7 74.4 

20 0-12 1010.5 167.5 73.9 
12-24 1010.0 168.3 74.0 

21 0-12 1011.4 166.8 74.1 
12-24 1014.2 74.1 

22 0-12 1017.6 166.3 73.9 
12-24 1022.4 165.5 74.0 

23 0-12 1026.8 167.6 74.5 
1028.0 166.4 74.8 
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cosmic ray 
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Lat. Long. 
50.3°N. 58.3°E. 
50.3 59.1 
51.6 62.3 
53.7 65.0 
55.7 68.5 
57.7 71.2 
60.4 72.4 
63.7 75.1 
65.1 76.2 
67.8 78.9 
70.4 81.4 
73.1 82.6 
75.9 84.9 
78.0 84.9 
80.0 75.0 
81.8 59.4 
82.9 47.8 
83.1 44.5 
83.1 44.5 
83.1 44.7 
81.8 
81.5 63.4 
81.3 63.3 
80.7 63.6 
81.2 63.8 
81.4 65.6 
81.2 69.3 
81.7 74.5 
82.3 73.4 

2.0 72.5 
82.1 79.8 
82.0 87.6 

2.9 97.7 
83.2 98.1 
83.2 98.2 
83.2 98.2 
83.2 98.2 
83.2 98.2 
83.2 98.2 
83.2 98.2 
83.1 97.7 
82.5 94.9 
82.0 85.4 
81.8 79.2 
81.2 76.3 
81.0 68.4 
82.2 53.6 
83.2 44.7 
83.2 44.7 
83.2 44.7 
83.2 45.1 
82.9 47.2 

2.6 42.6 
83.0 47.6 
83.1 43.7 
83.1 42.9 
82.7 39.3 
82.8 39.6 
83.0 38.4 
83.0 37.4 
82.9 35.1 
82.8 39.8 
83.2 42.7 
83.8 36.9 
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(A) Labrador (continued) 


Proposed 








| 
| 
| 


Corrected cosmic ray 
Date Time, Pressure, results, Lat. Long. Geomag. geomagnetic 
Gus. mb. counts ‘hour Lat. —_—_— 
+64 Lat Long. 
Aug. 24 0-12 1026.5 166.7 74.9°N. 101.9°W. 82.1°N. 84.2°N. 30.3°E. 
12-24 1021.3 168.4 74.9 104.2 81.5 84.4 24.5 
25 0-12 1014.7 167.7 74.9 105.7 81.3 84.5 20.7 
12-24 1009 .7 167.3 74.9 108.3 80.7 84.6 13.7 
26 0-12 1012.9 167.1 75.0 108.7 80.7 84.7 $2.4 
12-24 1017.5 165.5 74.9 108.8 80.7 84.6 12.3 
27 0-12 1017.4 166.1 74.7 109.3 80.4 84.4 10.7 
12-24 1012.9 167.6 74.5 110.0 80.1 84.3 8.6 
28 0-12 1006.5 167.1 74.3 110.7 79.9 84.1 6.6 
| 12-24 1004.2 165.4 73.9 110.9 79.6 83.7 5. 8° 
| 29 0-12 1000.6 167.5 73.5 111.4 79.0 83.3 4.4°E. 
12-24 1000.9 166.1 73.1 112.3 78.5 82.9 2.1. es 
30 0-12 1000.2 167.8 73.2 113.6 78.4 83.0 0.9°W., 
12-24 1003 . 4 167.4 73.3 115.3 78.1 83.1 5.0°W. 
31 0-12 1006.8 168.5 73.2 116.2 78.0 83.0 7.1°W 
2-24 1008.4 167.9 73.3 115.5 78.2 83.1 5.5 
Sept. 1 0-12 1012.1 167.7 73.4 115.6 78.2 83.2 5.8 
12-24 1014.2 169.2 73.4 114.8 78.3 83.2 3.9 
2 0-12 1016.6 167.6 73.3 115.7 78.1 83.1 6.0 
12-24 1018.0 168.2 73.3 116.3 Tae 83.1 7.4 
3 0-12 1019.7 167.3 73.3 116.4 77.8 83.1 7.6 
12-24 1020.2 168.5 73.2 116.4 C66 83.0 7.6 
4 0-12 1021.0 169.3 73.2 117.0 77.6 82.9 9.0 
12-24 1019.1 167.7 73.0 117.3 77.4 82.7 9.5 
| 5 0-12 1012.8 169.1 72.9 117.4 77.4 82.6 9.7 
| 2-24 1008.5 168.3 ta.¢ 117.8 77.2 82.4 10.4 
| 6 0-12 1009 .3 167.6 72.4 118.8 76.8 82.1 12.4 
12-24 1011.6 167.4 71.6 120.1 75.9 81.2 14.4 
7 0-12 1013.5 166.5 71.4 121.7 75.4 81.0 7.4 
12-24 1015.7 167.6 71.4 121.7 75.4 81.0 17.5 
8 0-12 1017.7 167.7 71.4 121.7 75.4 81.0 17.5 
12-24 1018.5 168.2 71.4 121.7 75.4 81.0 17.5 
9 0-12 1017.3 169.2 71.4 121.7 75.4 81.0 17.5 
12-24 1015.0 169.0 71.4 121.7 75.4 81.0 17.5 
| 10 0-12 1013.8 169.0 71.4 121.7 75.4 81.0 17.5 
| 2-24 1010.0 168.6 43.2 122.4 75.2 80.7 18.7 
| 11 0-12 1005.3 169.0 71.8 126.9 74.8 81.0 28.3 
12-24 1011.7 169.5 71.8 126.3 74.8 81.1 27.1 
12 0-12 1016.5 169.3 71.7 127.0 74.6 80.9 28.3 
12-24 1017.7 167.7 71.6 130.0 74.0 80.6 33.8 
13 0-12 1015.2 168.6 71.4 133.5 73.1 80.0 39.7 
2-24 1022.0 167.6 70.9 137.6 42.2 79 415.5 
14 0-12 1025.0 166.9 70.1 139.2 71.0 78.1 16.5 
2-24 1018.4 164.2 70.2 141.2 70.8 78.0 19.8 
15 0-12 1015.8 165.4 70.4 144.4 70.3 74.8 54.8 
12-24 1015.0 166.8 71.3 150.6 69.7 77.6 65.2 
16 0-12 1008.3 166.5 73.0 155.0 70.4 78.4 75.1 
12-24 1002.8 168.4 74.1 158.4 70.8 78.7 82.2 
17 0-12 1001.7 167.9 74.5 159.1 71.0 78.9 84.3 
12-24 1000.7 167.7 73.8 158.2 70.5 78.5 81.1 
18 0-12 1002.7 166.3 74.5 158.1 71.2 79.0 83.2 
12-24 1013.9 165.9 74.6 160.1 70.9 78.8 85.7 
19 0-12 1019.2 166.3 74.2 160.7 70.4 78.4 85.0 
12-24 1020.6 167.6 73.8 160.6 70.6 78.1 83.7 
20 0-12 1019.4 167.2 72.3 161.9 68.6 76.6 81.1 
12-24 1016.9 167.6 70.4 165.4 66.3 74.4 80.9 
21 0-12 1012.0 168.7 68.0 167.7 64.0 71.9 79.5 
12-24 1007.5 167.4 65.7 168.3 61.7 69.7 77.2 
22 0-12 1002.4 167.3 63.5 168.1 59.6 67.7 74.6 
12-24 998.4 168.7 60.7 168.1 57 65.1 72.2 
23 0-12 994.4 168.5 57.5 167.0 54.3 62.3 68.7 
12-24 993.0 165.3 54.7 165.3 52.0 59.8 65.2 
24 0-12 993.5 164.6 53.6 160.8 51.9 59.5 59.6 
12-24 996.9 166.2 53.0 155.9 52.3 59.6 53.7 
25 0-12 1002.4 164.1 52.5 150.7 52.9 59.8 17.4 
2-24 1013.5 167.6 51.7 145.8 52.9 59.6 41.2 
26 0-12 1022.7 163.0 51.1 141.7 53.2 59.4 36.1 
12-24 1027.0 164.5 50.4 136.1 53.7 59.2 29.0 
i 27 0-12 1026.2 164.4 49.5 131.3 54.0 58.7 22.8 
| 12-24 1019.7 164.7 48.7 126.7 54.0 58.2 17.0 
} 28 0-12 1017.2 165.0 48.3 123.5 54.1 57.9 12.9 
| 12-24 1019.6 165.2 48.3 123.3 54.3 57.9 12.7 
20 0-12 1019.1 163.5 48.3 123.3 54.3 57.9 12.7 
| 12-24 1019.3 165.4 48.3 123.3 54.3 57.9 12.7 
| i eg 
| 
I 
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(A) Labrador (concluded) 





























Proposed 
Corrected cosmic ray 
Date Time Pressure, results, Lat. Long. Geomag. geomagnetic 
U.T mb. counts/hour Lat. —_——_——_ —_ 
+64 Lat. Long. 
Oct.. 1 0-12 1001.9 163.9 48 .3°N. 54.3°N. 9°N. 12.7°W. 
12-24 1007 .9 162.9 48.3 54.3 9 12.7 
12 0-12 1011.0 163.4 48.3 54.3 9 33.2 
12-24 1013.8 164.2 19.2 54.9 8 12.5 
13 0-12 1017.9 161.8 49.2 55.3 8 12.5 
12-24 1024.7 163.5 49.2 55.3 8 12.5 
15 12-24 1014.3 164.8 49.2 55.3 58.8 12.5 
16 0-12 1010.3 163.5 49.2 55.3 58.8 12.5 
12-24 1007 .0 163.4 48.8 54.8 58.4 12.6 
17 0-12 1009.5 162.2 48.1 53.7 57.7 14.3 
2-24 1016.9 161.2 46.4 52.0 55.9 14.6 
18 0-12 1016.9 6 44.0 49.7 53.5 14.2 
12-24 1015.6 0 41.2 46.9 50.8 14 
19 0-12 1016.1 146.9 38.4 4 44.3 48.0 i3.1 
12-24 1021.7 146.4 37.7 om 44.0 47.3 11.6 
20 0-12 1020.6 145.7 37.8 5 44.1 47.4 10.9 
12-24 1020.5 146.8 37.8 5 44.2 47.4 10.9 
21 0-12 1019.0 146.3 37.8 5 44.2 47.4 10.9 
12-24 1019.2 147.8 37.8 5 44.2 47.4 10.9 
22 0-12 1019.1 147.0 7.8 5 44.2 47.4 10.9 
12-24 1020.1 147.1 37.5 5 44.1 47.1 10.9 
23 0-12 1017.5 141.9 35.2 121.8 41.9 44.9 9.9 
12-24 1015.0 136.4 32.7 120.2 39.5 42.4 8.0 
24 0-12 1014.8 130.4 30.1 118.6 37.2 39.8 6.1°W. 
12-24 1015.6 124.8 27.9 116.4 35.6 37 .7 3.6°W. 
25 0-12 1016.3 118.5 25.8 114.0 33.0 35.6 0.9°W. 
12-24 1016.0 112.7 23.8 111.8 32.1 33.6 pg : 
26 0-12 1013.6 108.8 21.8 109.6 29.5 31.6 3.9°F 
12-24 1013.1 105.9 20.0 107.5 29.0 29.7 6.2°E 
27 0-12 1012.1 106.1 18.7 104.8 28.0 28.4 9.0 
2-24 1012.5 106.8 17.4 102.2 oF. 27.0 11.8 
28 0-12 1011.9 102.9 14.9 101.9 24.6 24.5 12.0 
2-24 1011 100.1 12.8 101.4 22.6 22.4 12.4 
29 0-12 1012.3 98.4 11.6 98.5 21.6 2). 15.4 
12-24 1011.9 99.5 10.7 96.1 20.9 20.0 18.0 
30 0-12 1011.2 99.4 9.9 94.5 20.3 19.2 19.5 
12-24 1010.5 97.8 9.2 93.0 19.8 18.4 21.0 
31 0-12 1011.7 100.5 8.4 90.4 19.2 17.4 23.7 
12-24 1011.7 99.2 7.5 7.5 18.4 16.3 26.6 
Nov 1 0-12 1011.8 98.6 Tek 84.7 18.2 15.7 29.5 
12-24 1011.9 100.4 ak 82.0 18.4 15.4 32.2 
2 0-12 1011.2 101.1 7.5 80.1 18.7 15.7 34.2 
12-24 1011.6 103.5 8.8 79.5 20.1 16.9 35.0 
3 0-12 1011.1 103 .7 8.9 79.6 20.3 17.0 34.9 
6 0-12 1009.9 106.6 10.2 78.8 21.6 18.2 35.8 
12-24 1009.2 106.7 10.8 then 22.2 18.6 37.6 
7 0-12 1010.1 109.5 11.4 75.3 22.9 19.0 39.6 
12-24 1011.2 112.0 12.2 73.5 23.7 19.6 41.5 
8 0-12 1011.8 113.3 12.6 71.6 24.1 19.8 43.5 
12-24 1011.8 113.3 12.7 69.5 24.2 19.7 45.7 
9 0-12 1012.4 111.6 12.6 67.2 24.1 19.3 48.0 
12-24 1012.1 109.2 12.4 65.1 23.8 18.8 
10 0-12 1012.4 107.4 12.2 62.9 23.6 18.3 
12-24 1013.0 105.6 12.1 61.8 23.4 18.1 
11 0-12 1013.5 105.9 12.1 61.7 23.4 18.1 
12-24 1014.4 106.7 412.1 61.8 23.4 18.1 
12 0-12 1014.3 105.8 12.1 61.7 23.4 18.1 53.6 
12-24 1013.6 107.3 12.2 61.8 23.5 18.2 53.5 
13 0-12 1013.8 112.6 13.7 62.0 25.0 19.7 53.5 
12-24 1013.5 120.4 15.7 62.6 a1 21.8 53.2 
14 0-12 1013.6 121.0 a7 .7 63.1 29.1 23.8 53.0 
12-24 1014.5 127.3 19.7 63.8 31.1 25.9 52.6 
15 0-12 1014.9 131.5 21.6 64.2 33.0 27.8 §2.5 
12-24 1016.8 136.5 23.5 64.6 34.9 29.7 52.4 
16 0-12 1018.3 142.4 25.3 64.8 36.7 31.5 52.5 
12-24 1018.5 144.2 27.3 64.9 38.7 33.5 52.7 
17 0-12 1019.2 149.7 29.2 64.8 40.6 35.4 §3.2 
12-24 1019.7 153.1 31.2 64.7 42.6 37.4 53.7 
18 0-12 1019.5 155.9 33.2 64.5 44.6 39.3 54.3 
12-24 1021.0 157.7 35.4 64.4 46.8 41.5 54.9 
19 0-12 1021.9 161.3 37.5 64.4 48.9 43.5 55.4 
12-24 1020.5 161.6 39.5 64.1 50.9 45.5 56.3 
20 0-12 1018.3 161.2 41.1 63.8 52.5 47.0 57.0 
12-24 1015.6 163.9 42.7 64.0 54.1 48.6 57.3 
21 0-12 1012.4 162.3 44.0 63.5 55.5 49.8 58.2 


| 
| 
| 
| 








(B) Atka 


ROSE ET AL.: LATITUDE EFFECT 





Date 


Jan. 
(1955) 





- w 


ao 


1 


~] 


ao - & 





12-24 


0-12 
12-24 
0-12 
12-24 
0-12 
12-24 
0-12 
12-24 


0-12 
12-24 
0-12 
12-24 
0-12 
12-24 


0-12 
12-24 


0-12 


Pressure, 
mb. 


1016. 
1013. 
1011. 


1011 


1015. 
1017. 


1018. 
1016. 


1014 


1014. 
1014. 


1014 
1013 


1012. 


1012 


1012. 


1012 


1013. 


1013 


1013. 
1013. 


1012.6 
1013. 


1014 


1013. 
1013. 
1014. 
1015. 
1016. 


1016 
1018 


1018. 
1019. 
1019. 
1018. 
1017. 
1015. 


1015 


1016. 
1018. 
1021. 


1022 
1024 
1024 
1024 


1022 


ROOMS 


Ne ND wo He Die 00 


Wh RDOWNOH Watt hoon 


wWOURaH 


COMUNMON ROS 


ANUS 


PN Who to 


cavcoue 


Corrected 
results, 
counts/hour 


+64 


167. 
4 


164 


166. 
.6 
160. 
153. 
148. 


161 


145. 
140. 
133. 
127. 
124. 
116. 
108. 
105. 


106. 
106. 


151. 
152. 


ONO 


mH OmONmn 


MANOR NERORMONWD NBOHUNDH CHOWNRBNE 


0909 NICO OO Ore 


7 
0 


eK De OOuWo Oto 


oe 


Lat. 


~WBOWNONIN WIWON 


ONce wo 


ARAMAWDANS WORmWOIo 


WewwWiwwwe CMOREYUN NROWRIERa 


.5°N, 


981 





Long. 


70.5°W. 
69. 
69. 
70. 
‘ean 
72. 
73. 


144. 


147 


150. 
154. 


160. 
164. 


WNABQA AW NROUWWWO WANAAH 


Gmenene ammente 


ROMS COOH WON 


168.5 


172. 


176. 
179. 


177.6 
177. 
175. 


174.8 


174 


174. 
174. 
174. 
174. 
174. 
174. 
174. 
174. 


@® G0 G0 G0 Go GO G0 Go GO 


Geomag. 
L 


at. 


> on 
> res 
No WEIN WO 
yd 


& 


to 
-~ 
tor Cr be SIO he OD 


oa 
CSCO Uwonn 


t 
ANAWOwWAN CH WHODD 


wo 
to 
2 


> 
FADARAG HG 
WWOWWHWWH NWOORD COmMMmUwarIn 


Proposed 
cosmic ray 
geomagnetic 
Lat. Long. 
49.3°N 50.0°E. 
48.0 50.9 
45.4 50.0 
42.9 48.4 
40.4 46.6 
37.9 45.3 
35.5 44.0 
32.6 42.8 
30.2 42.1 
27.5 41.7 
25.1 40.4 
23.1 38.8 
20.9 37.0 
18.9 35.1 
17.3 34.8 
17.0 34.9 
17.0 34.9 
17.0 34.9 
17.0 34.9 
16.8 34.8 
15.2 34.1 
13.6 31.8 
11.4 29.3 
9.8 26.9 
8.0 24.3 
6.4 22.0 
4.3°N. 19.5 
ao. See 
1.0°N 14.9 
0.7°S. 12.6 
2.6°S. 10.1 
4.5°S. 7.5 
6.1 &.1°E. 
7.9 2.6°E. 
9.6 0.2°E. 
11.3 2.3°W. 
13.2 4.7°W. 
14.8 7.4°W. 
16.2 9.8 
17.7 12.7 
19.1 15.1 
20.5 17.6 
22.1 20.1 
23.5 22.7 
24.8 25.3 
25.9 28.0 
27.3 30.8 
28.4 33.8 
29.6 36.8 
30.6 40.2 
31.3 42.4 
32.4 45.6 
33.4 49.6 
34.2 52.7 
34.8 56.3 
35.3 59.6 
35.5 62.1 
35.6 62.3 
35.8 62.5 
37.2 64.1 
37.6 64.5 
37.6 64.5 
37.6 64.5 
37 .6 64.5 
37 .6 64.5 
37.6 64.5 
37.6 64.5 
37.6 64.5 
37.6 64.5 
37.6 64.5 





982 CANADIAN JOURNAL OF PHYSICS. VOL. 34, 1956 





(B) A tka (contin ued) 





Proposed 
Corrected cosmic ray 
Date i Pressure, _ results, Lat. Long. Geomag. geomagnetic 
J mb. counts/hour Lat.  —-————————— 
+64 Lat. Long. 





wo 
° 
wn 


Jan. 5 : 1016. 
1017. 

1013. 

1018. 

1017. 

1012. 

‘ 1009. 

12-24 1006. 


0-12 
12-24 


159.2 .6°S. 64.5°W. 


64. 


t 


mim RON OOO Ir 
ORR Oe 


SNe CH MORAN 


COM wownNwor~ri< 


ho NOD 
OOmaAS 


oO w 
OR OOO RENATO 


© 


mon 
PENOOONEN CHAR AIS: 


toe toed 


WhO RP ROO NPR OMNOe we Owe aleKg 


CUNDNDOOSO CwOR YEN 


Om: D 


> yy G9 iia at 


N& 
SHEN END Ou 


DOMWMANS NDODUNWO 


Mm Or aN Re Cb ROO ROO tN 


wwor 


~nrwmowoco~- 
MONNAWHO RWNHOOHRE CUWRR2wWwh 


ho 0000000 N00 


ihe ae a eee 
wSOmMNIN 


SA NOAA hm 


ae 


COCRHROR NON E ERE 
Sto 


weon 


t 
nhs AO = 


a 
9 
2 
3 
4 
5 
6 
0 
2 
0 
5 
2.4 
ae 
4 
9 
8 
5 
8 
8 
8 
7 
a 
8 
5 
4 
1 
Be | 
2 
od 
5.1 
2.0 
2 


ao 
MR RROHKODOAO BOWS SCOMODUON 


PNOOCOHOR AOWDBNUABDRHwWHON COBD WNWOOAHL 
NH wo Re NOB ANROWNN when WONDOREE 


OH We RONN WH We DOM WERUINO 


WOOCRHEWNN HOBBS 


CHBNBOOH 
‘ CenwoHrao 





ROSE ET AL.: LATITUDE EFFECT 


(B) Atka (continued) 


Proposed 
Corrected cosmic ray 
Date Time, Pressure, _ results, \. Long. Geomag. geomagnetic 
we mb. counts/hour Lat.  -————_—_—_————— 
+64 Lat. Long. 


3 
" 
% 

° 
m 


166.6 .4°S. 30.0°W. 59. 
165. 5! 

170. 
167. 
166. 
164. 
164. 
165. 


166. 


mONDwWWo 
" D 
2RWS Wedtomrary 


165. 
167. 
166. 
166. 
168. 
166. 


167. 
165. 
164. 
163. 
165. 
166. 
162. 


DWeNUSOP WO MONDO 
DOAROCOW NWOMOWL 


HKANNOAThR CHAGZOON CHONBWEtd 
CrwwNwro 
MN RAODD RWW DOW BRDDR ORD 


AM AOM WIN ROUNSTCAGN NUONNK 


MHONSmMHO wots 


BS 00 He 0 00 Or a 
NAROW Rm 


Coote moO tn 


POO POANOHR NRWR UWNORRNE 


~ 
oy t 
APAAARAW NP hehe RNOWH—wH NOWD CIRWHODRS 


ie 
.4 
2 
2 
5 
6 
.4 
.0 
3 
.2 
4 
3 
0 
0 
0 
3 
4 
5 
.2 
4 
0 
6 
8 
3 
i 
6 
8 
q 
3 
2 


AARAABO VZOCGCOCOCOOUROH OHOO 


CANDO 


ot 


toh Wewiwme UASOe 


KCORNSCBA S000 
CE ERD CHDUNCHR SHER PLERREA BOOCOCOOCOBONBA UNO NWDK WORD OCOUTOHN 


COM PROTON WWW WHwWiWWr SCHON NNNWRNRO OUR NOoUenr, 
yo o 3 D : 3 4 . > 
MWe CHYUASONCO SOOO COOCOON NRE RRR RR won CHwwr WwODUNUDE 


mie CObmMowWaun wor 
NDS Whe RWWOR AAAD 


oorun 
ono 





984 CANADIAN JOURNAL OF PHYSICS. VOL. 34, 1956 





(B) ua tka (concluded) 
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THE QUANTUM MECHANICAL CALCULATION OF THE 
SURFACE ENERGY OF SODIUM CHLORIDE— 
A FIRST APPROXIMATION! 


By F. vAN ZEGGEREN? AND G. C. BENSON 


ABSTRACT 


Léwdin’s theory for the cohesive energy of alkali halides has been used to 
calculate a first approximation to the surface energy of a {100} face of sodium 
chloride. The value found, 187 ergs/cm.?, differs from the experimental value 
determined by measuring heats of solution (276 ergs/cm.?) but is about 18% 
higher than the figure obtained from a corresponding classical model. For 
comparison the surface energies of several other alkali halides have been 
computed. 


INTRODUCTION 


Previous calculations of the surface energies of ionic crystals, such as the 
alkali halides, have been based in most cases on various adaptions of the 
classical Born—Mayer lattice model. Considerable complication arises when 
attempts are made to allow for changes in the state of polarization of the ions 
and in the lattice spacing near the surface of the crystal. Although no entirely 
consistent calculation of the surface energy of a crystal based on models of 
this type has been made, a reasonable estimate of the result of such a calcu- 
lation for the {100} face of sodium chloride is about 112 ergs/cm.? at 0° K. 
This is obtained by combining the value of 189 ergs/cm.? given by Shuttle- 
worth (1949) for an undistorted crystal with the distortion energy of 
77 ergs/cm.? which Benson, Schreiber, and Patterson (1956) have computed 
for a somewhat different model. In an experimental determination of the 
surface energy (or, more correctly, the surface enthalpy) of sodium chloride 
at 298° K., Benson, Schreiber, and van Zeggeren (1956) found a value of 276 
ergs/cm.? It seems unlikely on the basis of specific heat measurements 
(Patterson, Morrison, and Thompson 1955; Morrison and Patterson 1956) 
that the difference in temperature can account for the deviation between the 
theoretical and experimental results. 

Recently van der Hoff and Benson (1954) computed the surface energy of 
a {100} face of lithium fluoride using a quantum mechanical model and ob- 
tained a value considerably higher than that given by classical calculations. 
In view of the discrepancy noted above for sodium chloride a quantum me- 
chanical calculation for that crystal is of interest. Unfortunately some of the 
simplifying assumptions introduced in the treatment of lithium fluoride 
become less plausible for sodium chloride. In the case of lithium fluoride the 
positive ions were considered as point charges and the electronic charge distri- 
bution of the lithium ion was taken into account only indirectly through the 
orthogonalization of the orbitals about the fluorine nucleus. For sodium 


1Manuscript received June 11, 1956. : 
Contribution from the Division of Pure Chemistry, National Research Council, Ottawa, 


Canada. Issued as N.R.C. No. 4045. 
2National Research Laboratories Postdoctorate Fellow 1955-. 
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chloride a more explicit consideration of the electrons of the sodium ion and 
of the contribution of exchange terms is desirable. 

In this paper a quantum mechanical expression is derived for the surface 
energy of a {100} face of a NaCl-type crystal. Since distortion (displacement 
of the nuclei and polarization of the ions) in the surface region greatly com- 
plicates the mathematical treatment, these effects are not considered in the 
present calculation. Thus the formula derived gives only a first approximation 
to the surface energy. The model adopted is essentially the same as that used 
by Léwdin (1947 a, 6, 1948) for calculating the cohesive energies of crystals 
composed of atoms or ions with closed shell configurations. In the next section 
some of Léwdin’s notation and formulae are summarized and where necessary 
altered and extended with a view to their subsequent application. Then an 
expression is developed for the change in energy of an infinite NaCl-type 
crystal when it is separated along a {100} plane into two semi-infinite hemi- 
crystals. Finally the surface energies of sodium chloride and several other 
alkali halides are calculated and compared with the results obtained from 
classical models. 


FORMULAE AND NOTATION 


Léwdin (1947 a, 6, 1948) has shown that the cohesive energy of a system of 
ions with complete electron shells can be written in the form 


(1) E=ty Ew 


where E,,: is a term associated with the interaction between two ions g and g’. 
The sum is taken over all ions g and g’; the superscript (|g) on the summation 
indicates that terms of the form Ey, are omitted. In evaluating E it is assumed 
that the nuclei are at fixed sites and the effects of thermal vibrations are neg- 
lected. It is convenient to separate E,,, into four parts, 


(2) Eq’ = En t+Ege +E +E, 

which are called the Madelung-, correction-, exchange-, and S-energy terms. 
These are defined in such a way that they are separately symmetrized in g 
and g’. Lundqvist and Fréman (1950) have noted that this separation is 
somewhat arbitrary and that only the total energy is of real physical signifi- 


cance. 
The Madelung term for two ions g and g’ with valencies ¢, and e, and 


separated by a distance agg: is 
2 

(3) Ew = € €g€9'/ Aq’ 
where ¢ is the electronic charge. 

The correction term is written as 
(4) Ege = Iga Tee 
where explicit formulae for J“) and J?) are given in Léwdin’s thesis (1948). 
This energy represents a correction to the Madelung term arising from the 
fact that the ions are not ideal point charges at the lattice sites but overlapping 
charge clouds. 
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The exchange energy is a term of quantum mechanical origin 
re exch 
(5) Egg = Cog’ 


where C,,: is given by 


(6) Coe = a C(u|v). 


The superscript (g) indicates that the summation extends over all occupied 
orbitals y, of the ion g. The qty 


vi (1) y aes (1) ¥ (2) 5 
7) Culp) = — of f eC H ve) tel?) asa, 
is the exchange integral between two il one in the state y, associated 
with ion g and the other in the state y, associated with ion g’. 


The last term of equation (2) is 


(8) Evy 7a Bog +Dag 
and is called the S-energy term since it depends on the overlap integrals 
(9) S= fv.(an — i, 
The two terms on the right of equation (8) are 
(10) By = LB 
v 

and 
(11) Dag! = 2 e. 

¥ » 


Léwdin has shown that the first three terms of equation (2) can be considered 
as resulting from pure two-body forces. However this is not true for the 
S-energy term. Here both B,, and D,, depend not only on the ions g and g’ 
but also on the positions of all the other ions in the crystal. According to 
Léwdin 


(12) By, = BY +Bi +BY 
where B{? is a many-body potential given approximately by the formula 
(13) By = 2e" Sun Do Ar Qrusr 
h= 
and 
(14) A, = > ~ P,(cos 8,-"). 


Q’’ 


P (cos 6,--) is the Legendre polynomial of degree h and 6, is the angle between 
@,, and @,,. The sum extends over all ions g’”’ excluding g. Explicit formulae 
for the two-body potentials B“? and B‘ are to be found in Léwdin’s thesis. 

In computing the cohesive energy of an infinite NaCl-type crystal Léwdin 
replaced the sum in equation (14) by the first four terms. Actually because of 
the crystal symmetry the coefficients A;, A2, and A; were all zero and it could 
then be shown that 
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x (1) 2 
(15) nw = 2e Sypeacl os 


In the present calculation of the surface energy this expression for B“ will 
also be used, but the approximation is not as good as in the case of the infinite 
crystal since with the decrease in symmetry the higher coefficients A;, A2, and 
A; no longer vanish. 

Quantities BY”, BY, and B®) are given by equations analogous to (10). 
For convenience F,,- is defined as 


(16) Fa ee 2 ee. 
v & 

then 

(17) Biv. = 2e?FyyAo. 


In equation (11) the quantity D,, is given by 


Qo'9' 9 Ago"". 


. i oi 
(18) D,, = —e san = 0S yee 5 
and accordingly 


fi 2 ") Eg" €,"" 
(19) Dog = —e Fay Sat rinse 
o* sr Qgg’") 


Qg'g’’ 


DERIVATION OF AN EXPRESSION FOR THE SURFACE ENERGY 


Consider an infinite NaCl-type crystal with nearest neighbor distance a. 
Suppose that the crystal is separated into two hemi-crystals with {100} faces. 
As stated in the Introduction it will be assumed that there is no distortion in 
the surface regions. After division of the crystal the two hemi-crystals are 
widely separated and there is no interaction between them. The surface energy 
of a {100} face can be written as 


(20) C= * { [U]Jrinar = [U]rnteiar } 


where w is the area associated with a pair of surface ions. [U] initia and [U]pnai 
are the energies of the whole system before and after separation and are stated 
per ion pair in the surface planes. 

The last three terms of equation (2) decrease quite rapidly with increasing 
distance of separation a,,, and in calculating the cohesive energy Léwdin 
retains only terms arising from nearest neighbors. Adopting this same approxi- 
mation and using the symmetry of all the terms in equation (2) the energy of 
the system can be written 


(21) Rice .s =e oo ES +E +E) 
g 9 g Qo’ 


where in all cases g is summed only over the positive ions of the system. In 
the second term the superscript (mg) indicates that the summation over g’ is 
restricted to the nearest neighbors of g. 

In one of the hemi-crystals a set of positive ions, one chosen from each layer, 
will be indicated by gi, go,...gs..., where g; is situated in the 7th layer 
counting from the surface layer (2 = 1) on into the crystal. The designations. 
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g: will be applied to these ions both before and after the actual division, i.e. in 
the initial as well as in the final state. With this convention 


(22) “— E+ ) ae Ess + Es + Exe] 


where now the summation g is over ve set of positive ions g;, one in the first 
layer, another in the second layer, and so on. 
For a positive ion g; in the initial state, the following relations are valid: 


9 


(23) SEM = —S aw, 

(24) rege 6(I I +5), 

(25) Sow se = 6C,_, 
7 

(26) Fee ae \-26 “F,_ay-+B? Bet 
7 

where 


“2 dus ( _ gy er reer tty 


(27) au 
rai 2 
= 1.74756 
is the ordinary Madelung constant. The subscript + — on C, F and other terms 
in the above equations is used to emphasize the fact that they refer to C,,:, 
F,,’, etc., where g and g’ are nearest neighbor ions. These terms depend only 
on the distance a and are the same in the initial and final state. From equations 
(22) to (26) the energy in the initial state is 
2 


(28) [U]oa = fant 0( 1241246, +B +BY 
2e” 
%E pou) (31 


The sum > 1 arises from the infinite number of ions in the set g;, each of 
i=l 


which makes the same contribution to [Uinta 
Proceeding in the same way it can be shown that for a positive ion g; in the 
final state the equations analogous to (23) to (26) are 


(29) z ” Bor = —“faw+(—1)' S(i)}, 
(30) PE = (P+ 12) (6-a(i-1)), 


(31) UE = Cy-(6—-8((—-1)), 


(32) Ener «| -2! Fy—(aut+(-1)' Ai) + B+ B® | 


x (6—8(1-1)) +5 Fy_(-1)' [AG-1) 427 @) 
+S (i+1)—(SA(i-1)+- S(a)) 6(¢-1)], 
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where the delta function 6(i—1) is included to cover the case of the first layer 
where there are only five neighbors instead of six. The sums “can be evaluated 
from the expressions given by van der Hoff and Benson (1954). In general 
foru > 0 


ot (u) 


I] 


Dre Ste Hse (- 18 Phe + (ls tu)? 
li l3 


le oe 
scr Dp ei SS 
tiodd = leodd V (li +12") {1 +exp[— rr (i? +1,"))} 
A particular value required in later work is -“(1) = 0.066011. 


From equations (28) to (32) it follows that 


(34) [U]nnar = [Ula += 2 (-1)*'.S (i) 


(33) 


ll 


4 M4740, +BO+BY—2 : fesse 


+2 . rf oy -y"-ro-70 |. 
i=1 
Using equation (34) and the formula (van der Hoff and Benson 1954) 
(35) Do (-1)'S% (i) = B 


=s>7r dr — eee MCT 4 soins 
tiodd l2odd V (li +12") {1 +exp[— rr (li +22") }}? 
= 0.06524647 


in equation (20) with w put equal to 2a? leads to the expression 


2 2 
(36) o= 1 je 4—| 41240, 4B24B%-2 © Fy aw | 


2a la , 
2 
+2 - F.164-S ay 


for the surface energy. 


CALCULATIONS AND DISCUSSION 


Equation (36) is analogous to the expression 


2 2 

(37) Eeonyton pate = Mates) a8 1240, +52-439-2© Faw | 

used by Léwdin to calculate the cohesive energy per ion pair of an infinite 
crystal. For sodium chioride Léwdin (1948) has evaluated the exchange, 
overlap, and other terms in this equation, employing the Hartree-Fock wave 
functions for the free ions as a first approximation for the orbitals in the 
crystal, and has tabulated the results as a function of the nearest neighbor 
distance a. In using these data it is most convenient to introduce the symbol 
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(38) Eve = 192 4124+C,_+B3B? +B —2 - Fins 


for the repulsive contribution to the cohesive energy per ion pair. In terms of 
this the surface energy at 0° K. can be written 


2 2 

(39) e.= SAE Big F,_(6B@-—S(1)) ¢ 
where F,_ is the same as Fy in Léwdin’s notation. To be consistent with 
Léwdin’s work the interionic distance is taken to coincide with that separation 
of the ions which gives a minimum value of the cohesive energy. For sodium 
chloride this is 5.20 Ly (1 Ly = 0.5292 X 10-8 cm.) and the values of F,_ and 
Exp are respectively 0.010058 and 730.8 X 10-5 Ey (1 Eq = 4.3583 X10-" erg). 
The surface energy obtained by substituting these in equation (39) is ¢ = 
187.0 ergs/cm.? This is about 90 ergs/cm.? less than the experimental result at 
room temperature. As noted in the Introduction it is unlikely that the tempera- 
ture difference can give rise to a deviation of this magnitude. There is, of 
course, the possibility that the experimental value is high because of the 
effects of surface structure. On the other hand, the present calculation gives 
only a first approximation. The inclusion of refinements corresponding to van 
der Waals interactions, non-nearest neighbor repulsions, and distortion in the 
surface region might significantly alter the result. 

For comparison a value of the surface energy of sodium chloride was calcu- 
lated for a classical model corresponding to the present quantum mechanical 
one as closely as possible. This assumed no distortion, simple electrostatic 
forces, and an exponential repulsive potential b exp(—a/p) between nearest 
neighbors. The constant p = 0.345X10-° cm., and } was adjusted to give a 
minimum of the cohesive energy at 2.75 A. This calculation gave « = 159 
ergs/cm.” which is only about 30 ergs/cm.? less than the quantum mechanical 
result. This is quite interesting since the quantum mechanical model did not 
employ any direct experimental information about the sodium chloride crystal. 

The surface energies of four other alkali halides were calculated using 
equation (39). In the case of lithium fluoride it was necessary to compute the 
integrals S,, for the overlap of the 2p and 2s orbitals of the fluorine ion with 
the lithium 1s orbital. The results are summarized in Table I together with 
the values calculated from the classical model described above. 








TABLE I 
a Quantum value Classical value 
Crystal inA inergs/cm.?_ in ergs/cm.? 
NaCl 2.75 187 159 
LiCl 2.68 177 166 
KCl 3.09 184 128 
NaF 2.29 219 214 
LiF 2.40 145 195 


In the first four cases the quantum model yields a value somewhat higher 
than the classical treatment. This behavior is reversed for the last crystal 








992 CANADIAN JOURNAL OF PHYSICS. VOL. 34, 1956 


lithium fluoride. It is doubtful if any physical significance can be assigned 
to the present lithium fluoride calculation since the value of the interionic 
distance a = 2.40 A used is considerably more than that observed experi- 
mentally a = 2.00 A. Calculations for a = 2.00 A lead to a negative surface 
energy. 

The value of the surface energy of lithium fluoride calculated by van der Hoff 
and Benson is 557 ergs/cm.? Apart from the approximations already noted 
in the Introduction the model used by these workers differs from the present 
one in that possible covalent bond character and a simple form of distortion 
near the surface were considered. 
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NOTES 


A NOTE ON 
THE EFFECT OF THE RESIDUAL GAS PRESSURE UPON THE SPACING OF 
MASS SPECTROSCOPIC DOUBLETS* 


By NEIL R. IsENoR,f G. RONALD BAINBRIDGE, PuHiLip C. EAstTMAn,f 
AND Henry E. DucKWorRTH 


The mass spectroscopic comparison of atomic masses is usually accomplished 
by studying ‘“‘doublets’’. A doublet is a pair of mass spectral lines produced 
by two species of ion with specific charge almost, but not quite, equal. Thus, 
for example, a doubly-charged ion, which appears on a mass spectrum at a 
point corresponding to half its mass, will frequently form a doublet with a 
lighter singly-charged ion. If the mass of one of the doublet members is known, 
then the mass of the other may be computed from a knowledge of the doublet 
spacing and the dispersion of the mass spectroscope. 

It was demonstrated by Dempster and Shaw (1950) that a positive ion, in 
its passage through a mass spectroscope, may suffer significant energy loss 
without appreciably changing in direction. Furthermore, the extent of the 
energy loss was found to depend upon the nature of the ion. For example, in an 
original scale model of a Dempster (1935) double focusing mass spectrograph 
at a residual gas pressure of 10-> mm. Hg, 8000 volt NO* ions were retarded 
1.05 volts more than Nit* ions. This resulted in a reduction in the spacing of 
the NOt+-Nit+ doublet at mass 30, without serious deterioration in the sharp- 
ness of the two doublet lines. At higher pressures, when benzene was introduced 
to provide comparison lines, doubly- and triply-charged samarium ions were 
found to be retarded by as much as 16 volts more than were the corresponding 
hydrocarbon ions. Finally, also at the higher pressures, the lines of Sm+* were 
actually split. This was attributed to the existence of two different electronic 
configurations for the Sm*+‘ ions, one of which interacts much more strongly 
with the residual gas than does the other. While the line corresponding to one 
kind of Sm*4 ion remained as sharp as at the lower pressures, the line belonging 
to the other was badly broadened and indicated that these ions had lost as 
much as 58 volts more energy. Although this effect is potentially a serious 
source of error in the mass spectroscopic determination of atomic masses, it 
has never been the subject of any sustained investigation. To the authors’ 
knowledge, the only other published work on the subject, apart from that done 
by Dempster and Shaw, is a well-intentioned but rather inconclusive experi- 
ment by Duckworth, Johnson, Preston, and Woodcock (1950) in which the 
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spacing of the 4 Pt'®—Cu® doublet at mass 65 was studied in the pressure range 
2 to 10X10-*§ mm. Hg. In this case no effect was found. 

We have recently tried to ascertain if the 4 Xe!%-—Zn® and C3;H,O-Ni®8 
doublet spacings are sensitive to the residual gas pressure. The necessary 
doublets were photographed in a Dempster-type double focusing mass spectro- 
graph (Duckworth 1950), approximately 2.3 times larger than that employed 
by Dempster and Shaw. An accelerating voltage of 15,000 volts was employed. 
Nine good photographs of the 3} Xe'**—-Zn® doublet were obtained, two at 
3.5X10-*, four at 1010-*, two at 26X10-*, and one at 30X10 mm. Hg 
pressure. In the case of the C;HgO-—Ni** doublet, 13 good photographs were 
obtained, two at 5.5108, six at 14X10~-*, and five at 59X10-*° mm. Hg. 
The mass differences calculated from these doublets are shown in Fig. 1, in 
which the several values at each pressure have been averaged. The residual 
gas pressures were measured by means of an ionization gauge, calibrated for air. 

The data for the } Xe'?-Zn® doublet are shown in the upper half of Fig. 1. 
Within the precision of the experiment the doublet spacing does not appear to 
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Fic. 1. The }Xe%-Zn* and C3;H,O—-Ni® doublet spacings as a function of pressure. 
The numbers on the ordinate of the upper part of the figure should read 50.0, 50.5, 51.0, 
51.5, 52.0. 
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be pressure dependent, certainly not markedly so. The data for the C;H,O-—Ni°8 
doublet are shown in the lower half of Fig. 1. The evidence here presented 
suggests that the doublet spacing decreases with increasing residual gas 
pressure, the rate of change being approximately —0.035 m.m.u. per 10-° mm. 
Hg pressure. 

We realize that these experiments do not go very far toward the elucidation 
of this problem. However, they do provide some fragmentary confirmation of 
the fact that retardation of high velocity ions by gases depends upon the 
nature of the ions. Further, they suggest that not all doublet spacings are 
pressure dependent. Finally, they provide for the particular doublet 
C;H,O-Ni®® some quantitative evidence of the magnitude of the pressure 
dependence. 

It is planned to investigate this problem further with the aid of a large 
semicircular mass spectrometer (radius of curvature = 9 ft.), which has been 
recently completed in this laboratory. 
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